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PROJECT OVERVIEW 
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Advanced prop-fan propulsion systems which are being considered for 
passenger aircraft have noise characteristics quite different from turbo-jet 
systems they may replace. Rather than being primarily broadband much of the 
acoustic' energy of propellers is concentrated in tonals at the harmonics of 
the blade passing frequency (BPF) . Since the fundamental BPF is in the range 
of several hundred hertz, the tonal spectrum of the propeller induced noise 
is of concern from the viewpoint of passenger acceptance and comfort. 

The objective of this study is to examine analytically the structureborne 
contribution to the interior noise that is induced by the propeller wake acting 
on the wing. Analytical models have been developed to describe each aspect of 
this path including the excitation loads, the wing and fuselage structures, and 
the interior acoustic space. Because this is a feasibility study, the emphasis 
has been placed on examining a variety of parameters, and as a result different 
models have been developed to examine specific parameters. The excitation 
loading on the wing by the propeller wake is modeled by a distribution of 
rotating potential vortices whose strength is related to the thrust per blade. 

The response of the wing to this loading is examined using beam models. A model 
of a beam structurally connected to a cylindrical shell with an internal acoustic 
fluid has been developed to examine the coupling of energy from the wing to the 
interior space. The model of the acoustic space allows for arbitrary end 
conditions (e.g., rigid or vibrating end caps). Calculations are presented using 
these models to. compare with a laboratory test configuration as well as for 
parameters of a prop-fan aircraft. 

Several of the major findings of this study are summarized as follows: 

1. The unsteady-lift reaction of the wing to the propeller- wake field 

is a significant mechanism for transfer of vortex energy to structural 
vibration. 

2. In the absence of significant scattering from major structural 
discontinuities, the propagation loss of structureborne noise along the 
wing at the BPF is negligible and cabin noise levels are comparable to 
the source pressure levels for a cylindrical fuselage. 

3. Stringers that are small and in effect structurally integral with the 
skin have minimal impact on propagation along the wing. However, 
cut-off phenomena associated with a parallel array of discrete 
stringers could be created and exploited to enhance propagation loss. 

4. The interior acoustic levels on average are not appreciably altered 
by the presence of structural damping in the fuselage skin. 
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SUMMARY, CONCLUSIONS, AND RECOMMENDATIONS 


Advanced prop-fan propulsion systems which are being considered for 
passenger aircraft have noise characteristics quite different from turbo-jet 
systems they may replace. Rather than being primarily broadband much of the 
acoustic energy of propellers is concentrated in tonals at the harmonics of 
the blade passing frequency (BPF) . Since the fundamental BPF is in the range 
of several hundred hertz, the tonal spectrum of the propeller induced noise is 
of concern from the viewpoint of passenger acceptance and comfort. 

Propeller noise can propagate to the interior of the fuselage by paths 
that can be categorized as being primarily either airborne or structureborne . 
The airborne path is expected to be the dominant path, and it is estimated that 
the required insertion loss for this path provided by the fuselage must be 
approximately 55 dB at a BPF of 200 Hz. Measurements using conventional 
aircraft suggest that the contribution to the interior noise of the structure- 
borne path is not far below that of the airborne path. Consequently control of 
the structureborne path is likely to be important in achieving acceptable 
interior noise levels. 

The objective of this study is to examine analytically the structureborne 
contribution to the interior noise that is induced by the propeller wake acting 
on the wing. Analytical models have been developed to describe each aspect of 
this path including the excitation loads, the wing and fuselage structures, and 
the interior acoustic space. Because this is a feasibility study, the emphasis 
has been placed on examining a variety of parameters, and as a result different 
models have been developed to examine specific parameters. The excitation 
loading on the wing by the propeller wake is modeled by a distribution of 
rotating potential vortices whose strength is related to the thrust per blade. 
The response of the wing to this loading is examined using beam models. A model 
of a beam connected elastically to a cylindrical shell with an internal acoustic 
fluid has been developed to examine the coupling of energy from the wing to the 
interior space. The model of the acoustic space allows for' arbitrary end 
conditions (e.g., rigid or vibrating end caps). Calculations are presented 
using these models to compare with a laboratory test configuration as well as 
for parameters of a prop-fan aircraft. 
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The principal conclusions resulting from the analytical modeling and 
numerical calculations are the following: 

1. The unsteady-lift reaction of the wing to the propeller- wake field is 
a significant mechanism for transfer of vortex energy to structural vibration* 

For a fixed value of design thrust, and for points on the wing where wake loading 
is maximum, harmonics of the wake-induced pressure vary inversely with both the 
square of the propeller radius and the reduced frequency of the tonal. Predicted 
and experimentally obtained source levels are in reasonable agreement. 

2. In the absence of significant scattering from discrete stiffeners, 
e.g., ribs, or from other major structural discontinuities, the propagation loss 
of structureborne noise along the wing at the BPF is negligible and cabin noise 
levels are comparable to the source pressure levels for a cylindrical fuselage. 

3. For a rib to enhance significantly the propagation loss along a wing 
it must constrain both the translational and rotational motion of the skin. The 
cr ^^ er ^ a f° r achieving this objective appear to be realizable with minimum weight 
penalty. In practice rib effectiveness will be limited by rib elasticity which 
reduces its scattering strength and allows more efficient propagation across and 
coupling into it, both effects leading to an increase in cabin noise, the latter 
via the rib-spar-fuselage path. 

4. Sound radiation (airborne) associated with structureborne noise 
scattering from structural discontinuities along the wing is negligible compared 
with the direct acoustic path. 

5. Stringers that are small and in effect structurally integral with the 
skin have minimal impact on propagation along the wing. However, cut-off 
phenomena associated with a parallel array of discrete stringers could be 
created and exploited to enhance propagation loss. The criterion for effective- 
ness is given as k^L < tt , where k^ is the wavenumber of flexural waves in the 
skin and L the interstringer spacing. Also, the (spanwise) length of the 
stringer array need only measure a flexural wavelength in the skin and therefore 
entails minimum weight penalty. 

6. Spar vibration levels directly excited by the unsteady lift are less 
than those inthe skin near the midchord. Also, in the absence of stiffeners, 
either ribs or stringers, coupling of structureborne noise from the skin into 
the spars is minor. The extent of coupling from the skin to spars in the 
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presence of ribs has not been determined, and in view of the previous 
conclusions this skin-rib-spar-fuselage path may be the dominant one for 
typical wing configurations. 

7. The interior acoustic levels on average are not appreciably altered 
by the presence of structural damping in the shell. 

8. It is feasible to develop an analytical model of the structureborne 
path of the shaker driven laboratory model. However the noise predictions of 
a model that includes finite shell dynamics with frames, beam-shell inter- 
actions and structural damping, but ignores fluid loading, acoustic absorption 
and shell stringers do not compare favorably with laboratory measurements. 

9. The feasibility of developing an equivalent model of the structure- 
borne path along an actual wing- fuselage structure is questionable in view of 
the implied computer costs and/or running times. (For example the surface 

area of an aluminum fuselage of 3.05 m (10 ft.) diameter, 12.2 m (40 ft.) length 

-3 3 

and 2.5x10 m (0.1 in.) skin thickness measures roughly 10 flexural wavelengths 

at a frequency of 200 Hz. Using only a one-quarter wavelength mesh size and 

3 degrees of freedom per node a finite element analysis of the fuselage skin 

4 

alone would entail 10 unknowns exclusive of the modelling of the wing or the 
interior acoustic space.) 

10. Asymptotic or "canonical" analytical models of specific structural 
elements and features characteristic of actual wing-fuselage structures, e.g., 
stiffened skins, can be used to identify propagation phenomena that are 
potentially useful for diagnostic and noise control purposes. However, here 
it is particularly important to use such models in conjunction with a measurement 
program, either full scale or model, for confirmation. 

Recommendations for further work fall into two categories, these being 
analytical modelling and experimental measurements. Further development of the 
analytical models examined in this report is suggested in the following areas: 

1. In modelling the propeller wake, include the effects of, (a) a viscous 
vortex core, (b) three-dimensional unsteady aerodynamics, and (c) the dynamics 
of vortex paths in the presence of solid surfaces. 

2. In modelling the wing, develop and implement models that trace the 
skin— rib-spar— fuselage path accounting . for both rib elasticity and the spar- 
fuselage frame interaction. 
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3. In modelling the fuselage structure include the effects of, (a) internal 
stiffeners, and (b) either acoustic absorption along the boundaries (see 4 below) 
or the loading of the internal fluid on the structure. 

4. In modelling the interior acoustic space implement the theory in terms 
of the eigenfunctions appropriate to large boundary absorption (e.g., Bessel 
functions of the appropriate complex eigenvalues) . 

Recommendations for experimental studies directed toward characterizing the 
structureborne noise path are as follows: 

5. Using the existing laboratory model as a diagnostic tool for analytical 
concepts and models, examine the dependence of the interior pressure levels on 

a number of parameters including spatial location, structural damping, acoustic 
absorption, and discrete structural stiffeners. 

6. Develop new laboratory models of both wing and fuselage structures 
based on the scaling laws implied by. existing analytical models. These analytical 
models may also be used to assess scaling law conflicts, tradeoffs and compromises. 

7. The extent to which cut-off phenomena predicted for a stringer-skin 
geometry enhances the propagation loss along existing wings should be determined 
experimentally. Also the feasibility of further exploiting this mechanism for 
noise control purposes should be studied by means of a combined analytical- 
laboratory measurement program. 

8. Using a full scale aircraft and with excitation by shaker or by 
propulsion system when possible, (a) measure the structural transfer functions 
along the wing skin and spars; and (b) measure the transfer functions of 
structural response to interior pressure. Structural intensity measurements 
may be useful to examining the power flow along specific components of the wing. 
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I. 


INTRODUCTION 


Advanced prop-fan propulsion systems promise to provide fuel efficient 
flight at high subsonic speeds. The propellers in these systems having 
multiple skewed blades represent a noise source with characteristics quite 
from the dominant noise sources of current turbo— fan engines. In 
particular, rather than being primarily broadband, much of the acoustic energy 
of prop-fan propulsion systems is concentrated in tonals at harmonics of the 
blade-passing frequency. Since the fundamental blade passing frequency (i.e., 
propeller rpm times number of blades) is in the range of several hundred Hertz, 
the tonal spectrum of the propeller is of concern from the viewpoint of passenger 
acceptance and comfort. 

The paths by which propeller noise propagates to the interior of the 
fuselage can be categorized as being primarily either airborne or structure- 
borne. We define the airborne path as the path of propeller induced excitation 
that acts directly through the air on the fuselage. This includes propeller 
radiated noise and the aerodynamic near field of a rotating propeller located in 
close proximity to the fuselage. The airborne path is expected to be the 
dominant path by which acoustic energy excites the aircraft interior. The 
airborne insertion loss required by a fuselage to reduce the propeller tonal 
at the blade passing frequency to an 80 dB (A) (re: 20 yPa) interior level is 

approximately 55 dB (i.e., assuming 135 dB(A) free field level at 200 Hz). A 
variefc y means achieve a large sidewall transmission loss have been 

examined in the literature. Additionally, analytical models have been developed 
to predict the interior noise level in response to acoustic excitation of the 
fuselage. A recent review of many of these models is presented in Ref. 1. 

In general the parameters governing the transmission of airborne acoustic 
energy to the aircraft interior are more clearly defined than those that control 
the transmission of structureborne energy which is the subject of this report. 

We define the structureborne path as that path by which propeller induced energy 
propagates to the fuselage via other structures, notably the wing. Although 
useful for the purpose of discussion, the distinction between airborne and 
structureborne paths is not rigorous since structureborne energy can radiate and 
become airborne energy and vice versa. The sources of structural excitation at 
propeller related tonals include the interaction of the lifting surfaces with the 
propeller wake and engine vibration. The unsteady wake of the propeller acts on 
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the skin of the wing and of any rear control surfaces in its path. This 
vibration excitation can propagate along the skin and the primary wing structure 
to excite the fuselage. Engine vibration is transmitted through the engine 
mounts directly to the primary wing structure. Both of these excitations can 
transmit energy to the fuselage which in turn radiates noise to the aircraft 
interior. 

Structureborne noise has been cited by several investigators as a 

potentially important consideration for achieving acceptable interior noise 
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levels in propeller driven aircraft. While the structureborne path, which 

in fact may be multiple paths, is not expected to be the dominant means by which 
acoustic energy reaches the aircraft interior, its contribution may not be far 
below that of the airborne path, especially in view of the large noise reduction 
required for the advanced tuboprop aircraft. For example, measurements of 
interior noise m^de on a conventional twin engine turboprop aircraft (deHaviland 
Twin Otter) with both wing and fuselage wrapped to reduce excitation by the 
propeller wake suggest that the structureborne noise is no more than 10 dB 
below that of the aircraft path. ^ For a single engine aircraft with a 
reciprocating engine, measurements reported in Ref. 6 indicate that both airborne 
and structureborne contributions to the interior noise are comparable in level. 

In a more recent study, the importance of structureborne noise in a single engine 
aircraft has been further documented. 7 Analytical modeling of the structureborne 
noise path is currently needed to support and further evaluate the results of 
these empirical studies. Although the fuselage and interior acoustic space have 
been modeled to examine parameters important to the airborne path, 1 the 
excitations used in these models have typically been a single acoustic wavenumber 
pressure field. 

The objective of the present study is to develop an analytical model 
(albeit simplified) of the structureborne path including excitation loads, wing 
and fuselage structures, and interior acoustic space. The excitation loading 
on the wing by the propeller wake is modeled by a distribution of rotating 
potential vortices whose strength is related to the thrust per blade. Beam 
and plate. models are used to represent the response of the wing to this loading. 

A model of a beam elastically built into a cylindrical shell with an internal 
acoustic fluid has been developed to examine the coupling of energy from the wing 
to the interior space. The model of the acoustic space allows for arbitrary 
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end conditions (e.g., rigid end caps). Calculations are presented using these 
models to compare with a laboratory test configuration as well as for parameters 
of a turboprop aircraft. 
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II. THEORETICAL PREDICTIONS OF WING LOADING INDUCED BY A PROPELLER WAKE 
A. Introduction 

The near wake of a propeller in uniform flow consists of a complex 
system of trailing vortex sheets, each originating at a blade trailing edge. 

Due to self-induced velocities, these sheets quickly become unstable and roll 
up into a system of mainly concentrated tip vortices in a helical arrangement. 8 
Neglecting for the moment the effect of self- and mutually-induced velocities 
after the roll up process is complete, the trajectory of each vortex filament 
or vortex segment may be assumed as determined solely by freestream convection 
(in the propeller far wake, downstream of the wing and therefore not of concern 
here, this assumption is not valid and vortex structure and positions become 
chaotic) . • 


In the near wake, therefore, as a propeller tip approaches the horizontal 
plane say from above, its rolled-up tip vortex is convected over the wing, and 
an instant later, under the wing. Thus the vortex appears to "cut through" 
and, as a result, causes locally time- and spanwise-varying angles of attack 
and accompanying lift. The unsteady lift so generated becomes the source for 
noise borne by the wing structure into the cabin. 

Here we develop a model for the wing loads induced by the propeller wake, 
and compare predicted levels to those measured by Miller et al. for a 0.61 m 
(2 ft) diameter propeller in subsonic flow. Finally, estimates are made of 
vortex-induced wing pressures for an actual aircraft presently considered by 
NASA a possible candidate for the new prop— fan technology. When compared to 
typical estimates of incident airborne noise, the results of this last calcu- 
lation bear out the need for study of the structural path as a potentially 
significant competing mechanism. 


B. Propeller Wake Model and Wing Aerodynamics 


Figure II-l shows the near-wake model. A system of N straight, 
infinite potential vortices in rotary motion causes a downwash w on the wing 
plane given by 


w(t,y) 


r 

2uR 


n y-cosmt + — 3 
k=l . . 2 . 2iT(k-l) 


1+y -2ycos[flt+ 


N 


] 


(II-l) 
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where t is time and y , the spanwise distance measured from the propeller 
hub point, has been normalized by the blade radius R; Q and N denote propeller 
rotational speed and blade number, respectively. T stands for the strength 
of each vortex and may be determined from the propeller design thrust value 
via a lifting-line analysis performed later in the; chapter. 

We define the normalized downwash ft(t,.y) by 
~ - w 

W " r/2irR (II-2) 

and show typical signatures in Fig. n-2, where values of fl , and N have been 

chosen to match those in Miller's experiment. The y = 1.01 position corresponds 

to the measurement position in Ref. 3. The figure confirms that as y = 1 is 

approached the predicted downwash becomes more and more impulsive and that 

the period of the relatively complicated w sum is 2 tt/ 2N, with an associated 

fundamental (in Hertz) equal to the blade-passage frequency M) / 2 tt. Based on 

this we write the vortex-induced wing lift coefficient C (t,y) in terms of 

L 

harmonics C (y) ; 

Lm 

C L (t,y) . Ee |j 0 C Lm (y)a i "® t j- (II- 3 ) 

The aerodynamic model for the wing assumes the latter to be infinite in 
span (no three-dimensional root or tip effects) and as a flat plate in cross 
section. The flow is assumed aerodynamically two dimensional with each wing 
reacting locally to the input downwash field (classical strip theory) . Also, 
for applications of interest it will be found later that the airfoil reduced 
frequency based on the fundamental «Nb/U, where b is wing semichord, is not 
small relative to unity, and so the model to calculate C in terms of downwash 

~ LJU 

harmonics w m should contain high-frequency unsteady effects and should not be 
given merely by ■ 

27TW /U 
xn 

/l-M 2 

i.e., the quasi-steady assumption where U and M stand for flight speed and 
flight Mach number, respectively. Moreover, since flight speeds will be in 
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the high subsonic range, compressibility effects other than that given by 
2 -i 

the (1-M ) Prandtl-Glauert factor should also be included. 


For an airfoil in subsonic flow passing through a nonuniform field 

. 9 . 

given in terms of gusts, Amiet has derived such an unsteadiness-compressibility 

correction to the quasi-steady lift coefficient (kept in curly brackets in 

Eq. II-4 below) . Applying it in its relatively simple high-frequency limit 

(his Eq. 20) , one then has 
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C. Harmonic Analysis of the Propeller Wake 

The analysis to obtain the wake harmonics w' from w as defined by 

m 

Eqs. II-l thru II-2, begins with an identity involving two arbitrary complex 
numbers t = , z = re^ a (e.g.. Ref. 10). 


T-Z 


* * 
T -Z 


1-r 

l+r 2 -2rcos (ct-f5) 


(II-5) 


where * denotes conjugate. 


For r < 1 we may now express the left hand side of Eq. II-5 as 

9° m 

1+2 Z, r cosm(a-B) 
m=l 


Since r, the magnitude of z, is by definition positive, this relationship 
is immediately applicable to the range 0 < y < 1. . Thus letting y=r. Sit = a, 

-B = - 2l - ^ - , and performing some algebra one finds that 


y-cos f2 1 + 




l+y 2 -2ycos £ftt + — IT - ^ ' — 


<8 m-1 L . . 2Tt(k-l)l 

= - E. y cosm f! t + — 1 

m=l l N J 


(II-6) 


By a similar calculation for r > 1 the left hand side of Eq. II-5 takes ' 
the form 


-I 1 + 2 7 £ 0 - Sm( _ a-e) V 

\ "* 1 r m / 
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and for 1 < y < 00 we find that 


y-cos Qt 

2n (k-1) 1 

+ N J 

2 ^ 

_ 2tt (k-1) 1 

1+y -2ycos 

Qt + N J 


cosm Q t + 


m=0 


2TT(k-l) 

N 


m+1 


(II-7) 


If the cosines on the right sides of Eqs. II-5 and II-7 are now expanded out, 

and the sum E is interchanged with Z, the following two factors appear 
k in 


N r 2TTm(k-l) "1 5 • r 2irm(k-l) ~\ 

<ii cos [— J ' kii Sln L~ J ! 


the first of which is zero unless m is a multiple of N, and then it equals N. 
The second sum is identically zero. These conclusions may be used to find, for 
example, that the right side of Eq. II-6 becomes 

“ Nm-1 

-N Z, y cosmNQt 
m=l 

The simple manipulations used to obtain the series expression of w for 
0 < y < 00 may be adapted to the negative ranges -1 < y < 0, -°° < y <-l. One 
writes 


1-r 


1-r 


2 2 
1+r -2rcos(a-B) 1+r -2 (-r) cos (a-B+fr) 

and then replaces y by -y > 0 and — ^ — by 27r ^ + it in the previous results . 

The final series representations are now given for w = w(t,-°° < y < *>) : 

1 


w = N Z_ cosmNftt 
m=0 


Nm+1 


(- 1 ) 


Nm+1 


(-y) 


Nm+1 


1 < y < 


-00 < y < —1 


(II-8a) 


w = N Z_ cosnmiNQt 
m=l 


Nm-1 


-y 


, , . Nm , , Nm-1 

(-1) (-y) 


0 < y < • 1 
-1 < y < 0 


(II-8b) 


so that w = N/y 
m 


Nm+1 


for 1 < y < 00 , etc. Several points of interest should 
be noted. First, a steady field exists for |y| > l(m=0) which is algebraically 
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invariant with respect to the number of blades N. Since the m=0 harmonic is 
odd in y (proportional to ±l/|y|), this steady component represents a net 
moment about the propeller axis. The authors wonder if it has ever been a 
factor in wing design for steady flight. 

Secondly, increasing N causes a more rapid spatial decay of each m > 1 
harmonic from the N value at y-1. For very large N the vortex spacing becomes 
infinitesimal and physically a steady limit should be approached. Mathem- 
atically, in fact, only the m=0 steady term for |y| > 1 remains. In this 

limit the solution is consistent with the fundamental theorem of electrostatics 
which states that the internal field of a surrounding conductor is zero and 
that the field outside may be interpreted as due to a single equivalent 
centrally located charge. ^ 

Third and lastly, the initial expression in Eq. II-l yields an infinite 
value of downwash w at y = ±1 at times when vortices pass through the wing 
plane (the denominator becomes (y-1) (y+1) ) . As y -t ±1 therefore, the time 
history of induced downwash becomes highly impulsive (recall Fig. II-2) and is 
given by a series with tonals of constant level N (transform of a delta function) . 

D. Calculation of Vortex Strength 

Figure II- 3a (from Ref. 8) shows a blade airfoil section in rotary 
motion. As always, lift is defined as the force normal to the relative free- 

/ 7~2 2 

stream / U +(ftr) . We model the radial blade load distribution as indicated 

in Fig. II-3b, i.e., triangular in shape. At the. tip, where the radial load 

drops abruptly from its maximum value L . to zero, a vortex of strength T 

tip * v 

equal to the blade tip bound vorticity T b must be "trailed away" according to 
Kelvin's circulation theorem. To obtain the strength of the bound vortex in 
terms of tip lift , we apply the two-dimensional lifting-line relation 

p/u^+(nR) 2 r b = L fci . ( n-9) 

As previously discussed, upon passage through the wing the free vortex 
induces on the wing surface a time- and spanwise-varying angle of attack 
a(t,y) ~ w/U from which the vortex induced unsteady airload may be calculated. 
With defined in Eq. II-4 the high-frequency compressible model yields the 
following relationship for the m harmonic of vortex— induced wing pressure 


II-5 


L (y) = | u 2 c <y> (11-10) 

vm 2 Lm 

Substituting for r fa from Eq. II-9 into Eq. II-l, and then using Eq. II-l 


with Eq. II-4 into Eq. 11-10, we obtain 


L (y) = 

vm 


(?) 


(-i) 

— ftNmb 


Lastly, the tip lift L fcip must be related to the blade thrust T fe . 
Neglecting drag forces, the differential thrust dT and lift dL are related 
8 

through 


dT = dLcosB = dL 


U 2 +((ir) 2 


so that the total blade thrust is given by 


T b 


I*® 


Substituting dL/dr = L tip * (r / R) for the assumed triangular loading, and 
evaluating the resulting integral, one finally obtains that 


V R 


( 4 )'- ( 4 )’ •““(¥) 


Putting the above into Eq. 11-11 and defining the total propeller thrust T as 

T N, the desired harmonic L (y) of the vortex-induced, lifting pressure 
b vm 

distribution is found: 


L (y) = ~~2 
™ NR 2 


w (y) ! 

m . ± — 

■W {/‘-(Sr)’- (4)' ■“•‘(f)} 


ftmNb 
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For a fixed value of design thrust T, L (y=l) (the maximum for each 

vm 

harmonic) varies as the inverse of the square of the radius R, and as the 
inverse of QNmb/U, the "tonal reduced frequency". No other N dependence 
exists since the N factor in w^ (Eq. II-8) cancels that other appearing in 
the T/(NR 2 ) group of Eq. 11-13. 

As indicated in the final result, L depends also on Mach number M 

vm 

and on advance ratio S3R/U. However, the possibility of propeller-wake airload 
control through these parameters appears limited due to restrictions of actual 
practice on high subsonic cruise speeds U (and so M) , and on blade tip speed 
C2R. For well-known aeroacoustic reasons the latter is preferably kept also 
in the high subsonic range. 


E. Application of Wake Model to Two Cases of Interest 


We may now apply Eq. 11-13 to predict levels of the lifting pressure 
distribution on the wing due to the propeller wake and compare them to those 
measured by Miller, 3 whose laboratory model is shown in Figure II-4a. Also, 
a prediction is made using actual aircraft parameters. 


In Ref. 3 the vane stands vertically and is aligned with a propeller 
radius, i.e., it intersects the propeller axis if continued upwards in the 
sketch. Thus, the radial position given by r/R — 1.01 in the experiment 
corresponds to the spanwise position given by y = 1.01 in our theoretical model. 


The parameters for the laboratory propeller are the following: 


Power Coefficient: 

Advance Radio: 

Number of Blades: 

Blade Radius 
Frees tream Mach Number: 
Propeller Efficiency: 
Vane Semi chord: 


c = 1.84 
P 

J = 3.06 (S2R/U = .97) 

N = 8 

R = 0.30 m (1 ft.) 

M = .6 

^ e ff = .8 (private communication) 
k = 0.15 m (.5 ft.) (assumed) 


From these, values are calculated for thrust coefficient and total 

propeller thrust of 2.1 N (.48 lb) and 963 N (215 lb), respectively. For the 

first (m=l) harmonic at the spanwise position y = 1.01 (so that w = N) 

m 
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Eq. 11-13 yields 20 log 1Q (L^ / p ref ) = 154 dB (Pref = 20 yPa) vs, the experi- 

m 

mental value of about 143 dB. The rate of decay for the tonals may also be 
computed from Eq, 11^13 and compared to measurement. Taking m=8, the level 
of the eighth harmonic will be down by 23 dB (2Q log^^Cl.Ol /8) ) from that 
of the first harmonic. 

Both of these results are shown in Fig. II-4b superimposed on measured 

values for comparison. Fairly good agreement is indicated, though predictions 

appear somewhat higher than measurement at the lower part of the spectrum. 

We put off momentarily discussion of possible reasons for discrepancy and 

briefly present a similar calculation for an actual aircraft. 

12 . 

Re veil, et al. have listed operating parameters for a typical business- 

size propeller drive aircraft. These are: 

Tip speed = 243m/s (800 ft/sec) 

Cruise speed = .8M 

Propeller diameter = 2.2 m (7.2 ft) 

Number of blades = 8 

Cruising altitude = 9210 m (30,000 ft) 

Using a value of 302 m/s (995 ft/sec) for the sound speed at 9210 m 
(30,000 ft) , one calculates a flight speed of 242 m/s (796 ft/sec) , and so 
HR/U z 1. A value of b = 0.91 m (3 ft) is assumed for the semichord at the 
engine spanwise position. From Eq. 11-13 then, at the spanwise wing points 
of maximum vortex loading one finds that for the first harmonic (f = 283 Hz) , 

20 log [L (y=±l)/P =20yPa] = 137 dB 

10 v _ ref 

m=l 

a number comparable in magnitude to those measured and predicted for the 
airborne path (Ref. 1, Fig. 5) . 


Summary of Results 


When comparing the theoretical results predicted here to experimental 
results, the fundamental question should be asked whether a significant portion 
of the measured pressure can be attributed to the unsteady-aerodynamic lift 
mechanism of the present model. The propeller wake effects measured in Ref. 3 
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were apparently interpreted there as locally high values of stagnation 
pressures induced by free vorticity. The original purpose of the vane was 
probably to provide a physical mount for sensing instruments, not to act as 
an aerodynamic mechanism through which additional vortex energy could be 
channeled into the structure. 

The mechanism in our theoretical model is an aerodynamic "reaction" 
mechanism. Pressures predicted on the wing plane are due to the aerodynamic 
presence of the wing; they would vanish if the actual lifting surface were 
removed because they are the impermeable wing's lifting reaction to an incident 
unsteady flow. This resulting force driving the vane or wing is balanced by 
an equal and opposite net fluid force (Newton's first law) which in the 
immediate vicinity of the wing appears smeared over in the form of an acoustic 
near field. Since the theoretical model uses a flat-plate airfoil, those 
predicted pressures are antisymmetric with respect to the flight plane, so 
that on the two sides of the vane they are equal in magnitude but opposite in 
sign. 

The fairly good agreement between theory and experiment suggests that 
the lift effect is significant, and that it can therefore serve as an 
efficient coupler of propeller wake energy to structureborne vibration. It 
should be pointed out, however, that as estimated here these levels probably 
represent conservative estimates for the mechanism. The model does not 
account for finite vortex core size, nor does it account for mutual three- 
dimensional effects of wing sections in the vortex-induced spanwise-varying 
downwash field. Both of these may be expected to lower levels somewhat due 
to their smoothing action. 

Another possibly important neglected effect is that on vortex position 
and motion by an "image" system inside the wing. The model here allowed the 
propeller wake to cut freely through the wing plane, when in fact potential 
theory would call for a more complicated trajectory due to interaction with 
the solid surface. 
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Fig. II. 1 - Wake model for a four-bladed propeller, indicating horizontal 
plane of wing and vortex-induced downwash. 
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Harmonic 


Fig. II. 4 - (a) Laboratory model in Ref. 3. 

(b) Comparison of theoretically predicted pressure harmonics (o) 
for y = 1.01 to experimental measurements (solid line); for 
N = 8, n = 640 rad/sec, r/R = 1.01. 
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III. STRUCTUREBORNE NOISE MODELS 

A. Structureborne Noise Propagation Along a Wing Structure 

1. Introduction 

In the following sections the characteristics of structureborne 
noise propagation along a wing-like structure are analyzed, Basic concepts 
and functional analyses rather than extensive numerical modelling are emphasized. 

First the relevant characteristics of the vibration source, as modelled in 
Chapter II, are summarized. Propagation along the (bare) skin is then analyzed 
followed by a discussion of the air/fuel space between skins. Scattering and 
radiation from ribs is examined as is the potentially beneficial impact of 
discrete stringers. Finally spar dynamics are considered followed by a summary 
of results. 

2. Source Characteristics 


As described previously the assumed source of structureborne noise 
on the wing is the unsteady lift associated with blade tip vortices. This 
driving force, as modelled, has the following general characteristics: 

a. At a given spanwise position (y) measured from the propeller 
hub (y=0) the unsteady pressures on the top and bottom skins of the wing are 
equal in magnitude and opposite in phase. 

b. The unsteady pressure (lift) is antisymmetric about y=0. 

c. For y > 0, the unsteady lift at the m^ harmonic of the 
blade passage frequency N varies as (y/R) a where R is the propeller blade 
radius. The quantity a is given by 


a = (mN-1) 0 < y/R < 1 

= - (mN+1) y/R > 1 


(III-l) 


To account for the above three factors, we assume a normal pressure 
acting on the wing in the form 


P(y;w) =P (u)) f (y/R) 

o 

where P (w) = NL (1) as defined in Eq. 11-13, 
o vm 


(III-2) 


f (y/R) 


c 


= . -sgn(y/R) | y/R | ^ |y/R| 1 1 

= sgn(y/R) |y/R| -(rnN+1 * | y/R ] > 1 


(HI-3) 
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and 


sgn(y/R) = S 


-1 

v 


I y/R | > i 
I y/R | < l 


(III-4) 


The maximum lift p Q (w) is given in Chapter II where it was shown that for 
parameters typical of small aircraft p Q (w) ~ 137 dB re 20 pPa. 

3. The Wing as a Flexural Waveguide 

Wing structures consist of top and bottom skins with spanwise 
stringers, structurally connected to each other by means of a chordwise array 
of ribs and two (front and rear) spanwise spars (Fig* III-l) . The unsteady 
lift acts mainly over the surface area of the skins and the resulting vibrations 
are transmitted to the fuselage along the parallel skin-stringer and spar 
structures. 

We first consider the flexural response of the skins to the unsteady lift 
ignoring chordwise variations in both source and structure, that is we treat 
the wing as a one-dimensional (infinite) waveguide. The acceleration of the 
wing, a(y;w), in response to a distributed harmonic pressure field, p(y s ?tA)) , is 
given by (Ref. 13, Chapt. 6) 

00 

a(y;cj) = (-itj 2 / 4EIk 2 ) / p(y ; w) { exp [ ik I y-y I ] +iexp [-k J y-y I ] }dy (III-5) 

ri.s i 1 s r's s 


where to is circular frequency and the relevant material and geometric properties 
of the wing are expressed in the parameters E (Young's modulus) , I (cross- 


sectional moment of inertia) and k f (flexural wavenumber) , with 
k f = (ya) 2 /EI) ^ (l+irj/4) 


(III-6) 


where p is the wing mass density per unit length and n a structural loss factor 
to account for energy dissipation. (In Eq. III-5 stringers may be modelled 
only to the extent that they modify, on average, the parameters I, m and n.) 

Non-dimensionalizing Eq. III-5 to the inertial response of the wing under 
the peak pressure. 
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A (y ; a)) = ya ( y ; co) /P Q (co) = (-i/4) 


f (y g /3) {exp[i|y-y g | l+lexpL-jy-yJ ] }dy g (HI-7) 


where y = k^y and 3 = k^R. Thus, as modelled, at a given harmonic of the blade 
passage frequency the normalized acceleration is a function of the distances 
from the propeller hub to the locations of the wing tip and the observer 
measured in terms of flexural wavelengths. 

To ‘this point we have n6t related the parameters of the mathematical model 

to specific wing properties. However, it is useful to keep in mind the range 

of parameters of interest. Material properties are assumed to be those of rolled 

aluminum and we take E = 6.9x10^ Pa (10x10** psi) and p « 2.7x10** kg/m** 

-4 2 4 

(2.53x10 lb. sec /in ). Of concern is the fundamental propeller blade passage 

frequency, which as before, and corresponding to N=8, is taken to be 283 Hz, 

and its harmonics, f = u)/2tt = 283 m for integer m, R = 1.1 m (3.6 ft.) is a 

m 

typical propeller blade radius. For a skin of (effective) thickness h(in) 

Eq. III-6 yields a flexural wavelength 


A f (in) = 2ir/k f (in -1 ) = 2irl0 2 [h (in) /f (Hz) ] * = 37[h(in)/N]* 


(III-8) 


(Note: A f (in) = .025 A f (m)) 

-3 

For example taking 48 Pa (1 psf) for the weight density of the skin, h = 1.7x10 m 
(0.07 in.) and A^ - .25 m (10 in.) at the assumed fundamental blade rate frequency 
of 283 Hz. 


For illustrative purposes Eq. III-7 has been evaluated for the above sample 

parameters. The results are plotted in Fig. III-2 as a function of the non- 

dimensionalized spanwise distance y. In these calculations r\ was taken to be 
-2 

5x10 which is moderate. Also shown in Fig. III-2 is a plot of Eq. III-3 
showing the spatial decay of the source pressure magnitude from its peak at 
y - k^R ~ 27. 

Although not evident in Fig. III-2, as y => 0, A(y;w) goes to zero in 

accord with the assumed asymmetry of the source about y=0. For values of 

y < k f R, wave propagation associated with the nonuniformity of the pressure 

field produces an interference pattern in the response curve. The on-average 

level is roughly 10 dB down from that which would obtain from a locally reacting 

- 2 

inertial response to the maximum lift pressure (0 dB) . For k^R < y < 10 the 
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acceleration levels begin to decay exponentially due to structural damping. 
Finally if the graph was extended well beyond y = 10 , the curves for A(y ; u) 
and f (y/3) would eventually converge indicating that the local acceleration is 
dominated by the inertial response of the waveguide to the magnitude of the 
local pressure. 

In the following sections some of the cross-sectional details of the 
wing are considered, namely the entrapped air or fuel between skins, rib 
scattering and radiation, stringers, and spars. 

4. Entrapped Air/Fuel Between Skins 

inertia relative to that of the wing structure can be 
dismissed a priori. In addition its compressibility may be dismissed primarily 
under the assumption of identical top and bottom skins since as discussed 
ear ^ er unsteady lifts on the skins are out of phase, that is there is no 
net compressive force on the wing cross-section as a whole. For example, 
modelling the skin-air layer-skin cross-section as an acoustic waveguide 
terminated by the inertia of the skins and assuming equal and out of phase 
pressures on the skins, (Fig. 3a) either skin acceleration is given by 

ya (y; oj) 

p “ [2-itany] / [ 2-itany (1+ia) ] (III-9) 

with y = wd/c, a = pc/tou and where d represents the depth of the air cavity and 
pc the density and sound speed of the entrapped fluid. For parameters of interest 
a « 1 and the compressibility of the air layer may be dismissed for all 
frequencies. 

On the other hand consider the situation that exists for those wing 
sections containing fuel. In this .instance a > 1 and y < 1 for the frequencies 
of interest* and thus the fuel is effectively incompressible. The inertia of 
the fuel is large relative to that of the plating** causing a reduced (locally 
reacting) response of the skin covering these fuel tanks but also scattering of 


* Assuming (p<=) fuel = .75 (pc) water and d = 0.15 m (6 in.), a - 3xl0 4 /f(Hz) 
-4 

and y ~ 8x10 f (Hz) . 

** For the assumed parameters (pd) fuel /('ph) skin - 4d(in) >> 1. 
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incoming vibrations propagating along the skin from the fuel tank boundaries. 
This latter phenomenon is discussed in detail in the following section in 
connection with ribs. 

5. Rib Scattering and Radiation 


Ribs present structural discontinuities to the skin-stringer 
structure, modifying its impedance and wave propagation properties but also 
causing sound radiation. To examine the effect of ribs on wave propagation 
we first model a single rib as a locally reacting structure with drive-point 
impedance Z^_ (Fig. III-3b) . We locate the rib at a nondimensionalized distance 
y from the propeller hub. The wing acceleration becomes 


ya(y ; co) 

P (to) 
o 


CD 



f (y > g (y #y ) dy+F (to) g (y,y ) 
s s s r r 


(III-10) 


where 

g(y ,y ) = (-i/4)exp[i|y-y I ]+iexp[- |y-y | ) 
s s s 


(III-ll) 


and F is the nondimensionalized rib reaction 
r 

F r (to) = F r (to)k f /P o (to) . 


By imposing the constraint condition that 
a(y r ;to) = (ito/Z r )F r (to) , 


or 


ya(y ;to) 

■ — -r — = F (to) /Z , 
P (to) r r 


we may solve for the rib reaction 
-ff(y s )g(y r ,y s )dy s 


F r (to) « 


z 1 +g(y ) 

r r J r 


(III-12) 


(III-13) 


(III-14) 


(III-15) 


where 

Z" 1 = -itoy/k f Z r . (III-16) 

The wing acceleration in the presence of the rib may now be computed by 
substitution of Eq^. III-15 into Eq. III-10. The result i-s shown in 
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= °°) located one 


Fig. III. 4 for an infinite translational impedance rib (Z^ 
propeller radius from the wing tip (y = 26 ~ 54) . It is observed (compare 
Figs. III-2 and III-4) that although a null is indeed created at the rib 
location, beyond a quarter flexural wavelength away the effect is minor. It 
follows that the presence of an array of such ribs would be ineffective in 
reducing the vibrations reaching the fuselage unless the rib array spacing 
was less than a quarter flexural wavelength (A f /4 - 0.062 m (2.5 in.), for 
the illustrative example being considered) . 

Now let us increase the complexity of our rib model by allowing it to 
constrain the rotation as well as the translation of the wing. In this case 
the rib, by virtue of its rotational impedance Z^_, creates a concentrated 
couple (M^) in addition to a transverse force (F^) . The’ acceleration field is 
now expressed by 


pa (y ; o)) 

P (u>) 
o 


/ f < y. 


) g (y »y s ) dy g +F r (<*>) g (y >Y r ) +i« r (u) 3g (y ,y r ) /3y 


(III-17) 


where 


M • (u) = M (w)k‘/P (w) 
r r r o 


By imposing the additional constraint that 


3a(y r ;co)/3y = (iw/Z^M (u>) 


or 


y3a(y ;u>)/3y 

— r = -M (o>) /Z 

p (oj) r • r 


(III-18) 


we obtain 


M (oj) 
r 


/ f (y s ) 3g(y r ry s ) /9y[sgn(y r -y g ) ]dy t 
z r ’" 1 -i3 2 g(y r ,Y r ) /3y 2 


(III-19) 


with 


-* -1 3 * 

(Z ) = -iwy/k-Z 

r f r 


(III-20) 


w 

By setting both Z^_ and Z^ equal to infinity we create an effectively 
clamped boundary at y = y^ thus blocking completely wave propagation across the 
rib. Accordingly the computed acceleration levels for y > y^ (-27.) are caused 
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solely by that portion of the pressure field extending beyond the rib. And 
it follows that for this case the excess propagation loss across the rib is 
primarily a function of (y/R) , the distance of the rib from the hub measured 
in terms of the propeller blade radius, and should be approximated by 
20 log|y r /R| ^ niN+1 ^ . For the parameters being considered, m=l, N=8 and y r /R = 2, 
this expression yields -54 dB. Exact numerical results are shown in Fig. III-5. 

Figs. III-4 and III-5 were computed assuming a rib with effectively 
infinite translational and rotational impedances relative to those of the 
skin. Nevertheless, given the large skin compliance, this is readily achieved. 

The criteria for accomplishing this end may be deduced from the denominator of 
Eqs. III-15 and III-19 and are given by 

|Z r |/2/2 oip/k f = y > 1 (Ill-21a) 

and 

|Z r |/2/2 uiy/k^ = y > 1 

For example, assuming an inertial rib Eqs. III-21 are satisfied (y = 3) for the 
example being considered by a rib with a cross-sectional area of 6.2xl0~ 4 m 2 
(1 in. 2 ) and rotational inertia of 9.8xl0~^ m 4 (2.5 in. 4 ). 

Having considered the influence of single ribs the question arises as to the 
effect of a ribbing array. For example it is well documented^ 4 that a periodically 
spaced array of ribs produces "pass" and "stop" frequency bands within which 
flexural waves in the plating either propagate freely, that is unattenuated in 
the absence of dissipation or, suffer an evanescent decay. However, given the 
thin and, in turn, low impedance wing skin, it has been shown that the criterion 
for achieving a very large impedance mismatch at each rib location for both 
translation and rotation is easily realizable (Eq. III-21) . Also, for the 
problem at hand and as is the case for a single rib, the effectiveness of an 
array will be bounded by the direct excitation, it being distributed over the 
surface of the plating. Thus for the assumed parameters of interest, array 
effects (stop bands) are not considered to be a viable means for enhancing 
propagation loss along a wing beyond that implied 'by the multiple application 
of the above single rib analysis. Nor are array effects (pass bands) thought 
to be of concern in terms of short circuiting single rib effectiveness. 
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At this point it is informative to examine the dependence of the above 
predictions on N, the number of propeller blades. For this purpose the previous 
example was recalculated but for N=4 rather than 8. The frequency was held 
constant at 283 Hz. Results are shown in Fig. III-6 for the unribbed and fully 
constrained (in both translation and rotation) cases. Also shown in Fig. III-6 
as a dotted line is the spatial variation in the source pressure magnitude which 
is now observed to decay more slowly. However for the unribbed case, the 
normalized acceleration levels themselves are relatively insensitive to the 
change in N, especially the decay in level beyond y = 8. This may be explained 
as follows. For large values of mN the effective source area is concentrated 
around y = ±1, the locations of the wing tips. However, for the assumed 
parameters, the spatial decay in the source is still slow enough so that its 
effective' correlation length exceeds the characteristic (flexural) wavelength 
in the skin. To obtain a measure of the net source strength we integrate the 
pressure over one quarter flexural wavelength from the (closest) wing tip. 

Using Eq. III-3 we obtain an effective moment acting on the waveguide of 

(1+X /4R) 

2 ' f - 2 

M(w) = P(oj)R / (y -l)f(y /B)dy = P(tj)(X./4) 

J s s s f 

(1-X f /4R) (III-22) 

X f /4R « 1 

which is invariant to N. The (far field) waveguide response to this line 
moment is given by 

| A (y ; tu) | ~ (tt/ 4) 2 exp[-(y-g) n/4] . (III-23) 

Eq. 23 has been plotted in Fig. IIX-6 where it is shown to compare favorably 
with the far field portion of the response curve computed using Eq. III-7. 

While N has only a minor effect on the computed acceleration levels for 
the homogeneous structure, this is not the case in the presence of a strong 
structural discontinuity. For example, for the effectively rigid rib considered 
earlier, the transmission loss across it is roughly 50 dB when N=8 (Fig. III-5) 
but less than 30 dB for N=4 (Fig. III-6) . This difference in rib effectiveness 
is directly attributable to the slower decay in the unsteady lift for N=4 and 
in turn a higher excitation level at and just beyond the rib location. 
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Although the subject of this work is the propagation of structureborne 
noise along the wing it is important to keep track of the acoustic radiation 
associated with these vibrations as this radiation impinging on the fuselage 
is also a potential source of cabin noise. This combined structureborne- 
airborne noise path can short circuit potential structureborne noise control 
measures. 


To estimate this effect in connection with rib scattering we assume that 
the structureborne path has been blocked at some distance y^_ from the propeller 
hub. That is, at spanwise distances y > y , the wing acceleration is zero. 

The pressure radiated along the wing surface may be approximated by (Ref. 13, 

P- 94) , 


P r (y;ai) = (pW/2n) 



{exp[ik(y-y ) ] / (y-y )}a(y ; u))dy 
s s s s 


(III-24) 


provided that a = pc/com << 1 and (y-y^) >> W, where W is the average wing chord. 

In Eq* III-24 k = oj/c is the acoustic wavenumber. Rather than evaluate Eq. III-24 
exactly for the acceleration fields computed earlier let us consider a simpler 
case* 


The far field grazing radiation from a rigid piston centered about x=y=0, 
of rectangular dimensions 2L , 2L , and vibrating in an infinite baffle with 

* y 

acceleration amplitude a is (Ref. 13, p. 94) 


| p r (y-y r ;a) > I = tp/2ir(y-y r )]4a o L x sin(kL ) /k 


(I II- 2 5) 


We let 2L^ w the wing chord. Also, since for the example under consideration 

we are below the critical frequency, that is k £ » k, we allow 2L = A /4 or 

f y f 

one quarter of a flexural wavelength, it being assumed that radiation from the 
remaining portion of the vibrating wing is minimized because of phase 
cancellation. This yields 


P (y;w) | ~ pWa /4k -y 


y >> y. 


O' 'f J - -r 

Normalizing this radiated pressure to the maximum source pressure 

| p _ (y;w) /P (to) I - pWa /4k.yP (oj) 
r o O I o 

= ( 4/l2 k/k f ) -1 (p/p w ) (c w /c) (W/y) (ya_/P_(a))) . 


(III-26) . 


(III-27) 


o o 
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For example, using the physical parameters of our example and a frequency of 
283 Hz (m=l) 


|P (y;io) (P /co) | - 4xlO -3 (W/y) (pa /P (to) ) 
• r o o o 


(III-28) 


Taking pa /P (oj) to be roughly unity (see Fig. III-5 with y~y ) we see that, 
o o r 

ignoring spherical spreading in the form of (W/y) the radiated pressure is down 
a considerable 48 dB from the maximum source pressure. Given the Chapter I 
estimate of P Q (w) = 137 dB re 20 yPa for a typical maximum source pressure, 
it is concluded that this radiated field is negligible compared with the direct 
airborne path. 


6. Discrete Stringers 

To this point stringers have been modelled as being integral with 
the wing skins and thus accounted for only to the extent that they modify skin 
parameters yielding an effective stiffness or density. However, depending upon 
geometric details this approach may not be adequate. For example if the 
stringer impedance is large relative to that of the skin and if the stringer 
spacing is comparable to or exceeds the wavelength of freely propagating 
flexural waves in the skin then the stringer must be modelled in a discrete 
fashion. It is shown below that this situation provides a mechanism for 
potentially enhancing the propagation loss of structureborne noise along the 
wing over the bare skin, or integral skin-stringer case. 


Consider the structure shown in Fig. 111^7. The skin which is still taken 
to be of infinite spanwise extent is now taken to be not just a lifting line as 
before, but rather a lifting surface of finite chord, simply supported along 
its sides. These simple supports represent the stringers. Assuming an 
excitation pressure of the form F(x,y) = R(x)5(y), the plating acceleration, 
which is now a function of x and y, must satisfy the partial differential 
equation 


2 

a +2a +a -k^a 
yyyy yyxx xxxx f 


-W 2 [12 (l-v 2 ) /Eh 3 ] f (x) 6 (y) 


(III-29) 


The solution is in the form of a series 
a(x,y;u) = I a n (y) sin(mrx/L) 


(III-30) 


Taking the Fourier transform of Eq. Ill- 2 9 with respect to y yields 
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(Ill— 31) 


4 2 2 4 4 - 

[(k +2kV+k )-k*]a 

y n y n f n 


with 


f = -oj 2 [12(1-v 2 ) /Eh 3 ](2/L) / f(x)sin(k x)dx 

n In 


and where k^ - nn/L and k^ is the wavenumber corresponding to the y- transform 
variable. Inverting Eq. III-31 and assuming a uniform chordwise pressure of 
unit amplitude, that is f(x) =1, 


pa(x,y;io) = (-ik / 2 ) * (2/mr)sin(k x)G (y) 

r n-i,j,b n n 


with 


G (y) - ^exp[i(l-K 2 ) ^y]+i(l+K 2 ) *exp[- (1 +k 2 ) *y] (III-3: 

11 ii n n n 

where k = k /k _ and as before, y = k^y. 
n n f f 

Eq. III-33 may now be used as an influence function governing the wing 
response to the distributed lift defined previously but now in the presence of 
discrete stringers , 


ya(x,y ; w)/P (to) = (-i/2) , _ Z _ 

o n=l,3,5. 


(2/niT) sin(k 


f (y /P)G (y ) dy 
s n s s 


Although detailed calculations using Eq. 34 could be performed it is 

sufficient, for present purposes, to examine the nature of the solution. Of 

principal interest is that Eq. III-33 implies a cutoff frequency below which 
2 

K n > 1 and therefore the acceleration decays exponentially with spanwise 
distance from the source. The lowest cutoff frequency, corresponding to n—1, 
is given by < 1 » 1 or k f L = tt. 

For example, at one-half of this frequency k f L = tt// 2, and the minimum 
exponential decay (n=l) given by Eq. III-33 is exp[- ( tt// 2) y/L] . Thus, the 
decay at a spanwise distance from the source equal to only one stringer spacing 
(y=L) yields -20 dB of attenuation. 

Unfortunately this rather impressive effectiveness can be short-circuited 
by the flanking path created by the rigid body motion of the stringers. This 
places a restriction on the minimum size of the stringers, but fortunately only 
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on a spanwise length comparable to the stringer spacing. For example * a 20 dB 
reduction in this rigid body motion requires a stringer to plating mass ratio 
of the order of ten. 

Typical values of L may be estimated as follows. The cutoff criterion 
k f L = it yields f(Hz) = ir 2 10 4 h (in) /L 2 (in 2 ) (= 2. 5xl0 3 h (m) /L 2 (m 2 ) ) for an aluminum 

skin. Thus, taking h = 1.7xl0~ 3 m (.07 in.) as before, a 0.10 m (4 in.) stringer 
spacing yields a cutoff frequency of 430 Hz, well beyond the assumed fundamental 
blade rate frequency of 283 Hz. 

7. Wing Sweep and Spar Dynamics 

In previous sections, with the exception of the discrete stringer 
analysis, wing dynamics were modelled one dimensionally, that dimension being 
the spanwise coordinate. This can be justified on the basis of the source model 
which predicts an average lift invariant to the chordwise dimension and the 
implied assumption of a wide wing with zero wing sweep. As a consequence, the 
structureborne noise propagation characteristics of the wing are dominated by 
those of the skins. In this section we consider the potential structureborne 
noise path along the spars, this path being parallel with and coupled to that 
along the skins. 

The spars themselves not being directly exposed to the unsteady lift are 
excited via the ribbed skin structure. Since the overall surface area of the 
skin greatly exceeds that which is in direct contact with the spars, it is 
assumed here that spar motion is the result of propagating waves in the skin 
being coupled into the spars. This may be justified as follows. 

16 

Taking the chordwise variation in the unsteady lift to be of the form 
X (x) = it -1 [ (b-x) /(b+x) ] ^ with the front spar located at x = -b, and where 
b = W/2, the wing semichord as defined in Chapter II, the integrated pressure 
acting on the skin over the chordwise distance x^ > x > -W/2 becomes 

x 

c 

/W) + ( 1- ( 2x /W) 2 )*] . 
c c 

If the effective force on the spar (per unit spanwise distance) is approximated 
as that acting on the adjacent skin extending out a chordwise distance of A f /4 
then the ratio of this force to the overall force acting on the wing (per unit 
spanwise distance) is 


P (* ) = 
c c 


/ 

-W/2 


-1 


X(x)dx = (W/ir) [it/2+sin (2x 
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A f /W « 1 


P c (-W/2+X f /4)/p c (W/2) = (2 /it) (X f /W) * « 1 

To examine the skin-spar phenomenon we analyze the situation shown in 
Fig. III-7. A straight crested flexural wave propagating in a semi-infinite 
plate (skin) is incident on a beam (spar) of infinite extent at an angle given 
by <j>. This angle may represent the sweep angle of the wing or another source 
of skew. The beam serves as the termination of the plating and is capable of 
sustaining both flexural and torsional waves. However we will simplify our 
analysis by ignoring coupling between torsional motion in the spar and in-plane 
or membrane motion in the skin. The analysis of this problem is straight- 
forwardand follows closely that given in Ref. 17. (In Ref. 17 an infinite 
rather than semi-infinite plate is analyzed.) A unit amplitude harmonically 
varying flexural wave is defined by 

w^ (x ,y) = exp(ik f ysin<{>) exp (ik f xcos<{>) (III-35) 

where is the incident angle on the beam measured from the x coordinate axis. 
Grazing incidence on the beam (spar) is given by cj) = tt /2 and normal incidence 
by <{) = 0. The reflected wave takes the form 

2 ^ 

w r (x,y) = exp(ik f ysin<J»)+Rexp(-ik f xcos<{i) +R N exp(-k f x(l+sin <j>) ) (III-36) 

The resulting translational motion of the spar (x=0) is given by 

|w g (y) | = |w(0,y) | = IvMO^+w^Ojy) | = 1+R+R N (III-37) 

and the (torsional) rotation of the spar is 

| 8w g (y) /9x| = k x l 1-R+iR N l • (III-38) 

After expressing the interaction shears and moments in terms of the boundary 
motion the reflection coefficients R and R^ may be determined from continuity 
of displacement and slope. The! result is the following set of simultaneous 

equations. / 

i 

[A]R = B j (III-39) 

with ; 



R 


2 2 

“|S+l+sin <{>-vsin $ 

j 

R = 

_v 

, B = 

i 

/ 2 2 
ivJ cos({)-(cos (})+v sin 



/ 

L \ 

t 


I 

\ 

y 
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and 


[A] 

with 


2 2 

S+icos<}> (1+sin <f>-vsin <f>) 
2 2 

iUcos<{>+ (cos (j>+vsin <j>) 


S = (A^k^/h) [ ( c b /c^) 4 sin%-l] 
U = (Jk^/h) [l-(c t /c p ) 2 sin 2 <|>] 


2 i 2 

S-v(l+sin <}>) sin 

2 \ 2 2 
U(l+sin <f>) -(1+sin <j>-vsin (J>) 


(III-40) 

(III-41) 


2 A 1 

and where 0 ^ = 0 (K/J) and c, = (c(r ) , w) are the velocity of torsional and 
t s b g b ^ 

flexural waves in the beam and c = (c(r ) oj) t is the velocity of flexural 

P 9 P 


In the above equations (r )__ and (r ) are the radii of 

9 b 9 P 


waves in the plating 

gyration of the beam and plating respectively, and it has been assumed that 

the sound speed in both the beam and plating is c. The term c represents the 

s 

shear velocity in the beam and K and J its cross sectional torsional constant 
and polar moment of inertia respectively. Also h = /l2(r^) is the plating 
thickness and A^ the cross sectional area, of the beam. It is noted from 
Eqs. III-39 that with the exception of Poisson* s ratio v# the dependence of 
the reflection coefficients on the physical parameters of the system is solely 
through the flexural parameter S and the torsional measure U. Once again our 
interest will be limited to an analysis of the nature of the solution dictated 
by the mathematical model. 


First we observe that for grazing incidence (<j> = tt/ 2) , Eqs. III-39 
reduce to 

3 S— /2 v 


Lv 



R - 


“S ■ 

J 

-V 


_ v _ 


(III-42) 


yielding R = -1 and = 0. Thus for this , (zero wing sweep) case we confirm 
that there is no coupling into the £>eam (spar) . Although it is not particularly 
relevant to the problem posed, it is of general interest to explore the other 
extreme, namely normal incidence for which <{> = 0. Here Eqs. 39 reduce to. 
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(III-43) 


S+i 

S 


' R ‘ 


S-i " 





= - 


_iU+l 

U-l. 


-V 


_ 1-iU _ 


And it can be shown that under this circumstance 

f=“|W s (y) | 0 S » 1 

4 S,U « 1 

2 S « 1 , U » 1 

Physically , large values of S and U imply large translational and rotational 
beam impedances respectively, relative to the inertial impedance of the plating. 

More relevant to the issue at hand, we now consider the solution to 
Eqs. III-39 for arbitrary values of $ but a variety of limiting cases for S 
and U. For S >> 1, Eqs. III-39 yield the following solution for the reflection 
coefficients 


“ R “ 

= D _1 

~ U (l+sin 2 <J>) t (l+(l-v) sin 2 ^] -l“ 



“ R N“ 


2 2 

__ - [iUcos<f>+cos <}>+vsin <{>] 1. 


2 2 

_ iUcos(j)- (cos <}>+vsin 


(III-44) 


with 


2 \ 

D = U [ (1-fsin (f)) -icos<J)]-2 . 


If in addition we assume U >> 1 then Eqs. III-44 give R+R^ ~ -1 and therefore 
from Eq. III-37, 


L » 

S,U=s>°° s 


(y) = o 


(III-45) 


On the other hand if we assume U « 1 then Eqs. III-44 yield 

L 

Ss.JTu- 0 V y) = vsin ♦ • 


(III-46) 


Similarly it can be shown from Eqs. III-39 that 
1 ! 2 2 

U=S> J“£ =0 W s (y) =j (2+sin <(.)/ (l+sin <{,) . (III-47) 

Thus, in the asymptotic} limit of large translational and rotational spar 

i 

impedances, there is ncj translational coupling into the beam. However, if 


either the translational or rotational coupling factors is negligibly small 
then the coupling is of the order unity. 

The quantity S is small if either A k /h « 1 or sin<J> ~ (c /c ) . 

b f p b 

However, since for typical wing structures A fa k f /h » 1 only the latter 
condition is of interest. Taking the spar and skin to be of the same material 
(c p /c b ) = [( r g) p / ( r g ) ^ > an ^ therefore there is strong coupling into the spar 
for incident angles 


sin 1 [ (r ) /(r ) ] * 
g p g b 


(III-48) 


Similarly, since typically Jk^/h » 1 there is strong coupling into the beam 

r i 

for incident angles that satisfy sinA ~ (c /c ) or 

P t 


<$> s sin 1 [ (tt// 3) (c/c )k *h /X ] (III-49) 

n s p c 

where the geometrical factor k = K/J and X^ = 2 ttc/u) is the wavelength of 
compressional waves in the plating. Eqs. III-48 and III-49 are expressions for 
the angles at which there is wavenumber matching between flexural and torsional 
waves in the beam and flexural waves in the plating. 

In order to evaluate the degree of coupling between wing skins and spars 
based on the above model, one first estimates the quantities S and U. To do 
so we hypothesize an aluminum spar, rectangular in shape, of thickness t = 2h 

*3 • S p 

and height H = 100 h where the skin thickness h_ = 1.7xl0 _ m (.07 in.), 
s p P 

Further we use c/c = 1.7 for aluminum k = 4(t/H)^ for a rectangular cross 

18 S 

section, and consider a frequency of 283 Hz. This yields <$> f - 6 Q and <J> ~ n° 

and thus strong coupling due to wavenumber matching occurs only near normal 
incidence. Also A^k^/h =8.8 and Jk^/h = 14.2. Therefore, if we somewhat 
arbitrarily assume a wing sweep of 30° and take this as our angle of incidence 
then for the above parameters S,U > 10. From these large values little coupling 
is expected between the (bare) skin motions and the front and rear spars. 

Although beyond the present scope, it should be noted that the above 
conclusion may not hold in the event that the effective skin thickness is 
significantly larger than assumed or if there are large rib motions coupling 
well into the spars. Also, in this discussion the translational motion of the 
spar was taken as a measure of coupling while in practice spar rotations are 
also relevant. 
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8. Summary of Results 


The more salient results from Sections III.A.l through III. A. 7 
are summarized below. 

a. For a rib to enhance significantly the propagation loss along 
a wing it must constrain both the translational and rotational motion of the 
skin. The criteria for achieving this objective appear to be realizable with 
minimum weight penalty. In practice rib effectiveness will be limited by rib 
elasticity which reduces its scattering strength and allows more efficient 
propagation across and coupling into it, both effects leading to an increase in 
cabin noise, the latter via the rib-spar-fuselage path. 

b. Sound radiation (airborne) associated with structureborne 
noise scattering from structural discontinuities along the wing is negligible 
compared with the direct acoustic path. 

c. Stringers that are small and in effect structurally integral 
with the skin have minimal impact on propagation along the wing. However, cut- 
off phenomena associated with a parallel array of discrete stringers could be 
created and exploited to enhance propagation loss. For an illustrative example 
the criterion for effectiveness is given as k^L < tt , where k^ is the (frequency 
dependent) wavenumber of flexural waves in the skin and L the .interstringer 
spacing. Also, the (spanwise) length of the stringer array need only measure a 
flexural wavelength in the skin and therefore entails minimum weight penalty. 

d. Spar vibration levels directly excited by the unsteady lift 
are less than those in the skin near the midchord. Also, in the absence of 
stiffeners, either ribs or stringers, or extreme wing sweep, coupling of 
structureborne noise from the skin into the spars is minor. The extent of 
coupling from the skin to spars in the presence of ribs has not been determined, 
and in view of the previous conclusions this skin-rib- spar- fuselage path may be 
the dominant one for typical wing configurations. 

B. Wing-Fuselage Interaction 
1. Introduction 

In this section a simplified analytical model of a fuselage that 
is structurally connected to a wing is presented. In an aircraft the frame of 
the wing, namely the front and rear spars, is structurally connected to that of 
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the fuselage via a geometrically complex structure, shielded somewhat from the 
passenger compartment of the cabin. Thus structureborne noise levels in the 
spars will be attenuated to some extent as they propagate through the fuselage 
and eventually radiate as cabin noise. For an aircraft the detailed modelling 
of the wing-fuselage structure is beyond the scope of this report. Nevertheless, 
to examine some of the characteristics of wing-fuselage interactions and to 
simulate a laboratory experiment recently performed at NASA/Langley a simplifed 
analyticalmodel is studied. In this model the fuselage is represented by a 
uniform cylindrical shell while the wing is modelled as a freely supported beam. 
Coupling between the beam and shell occurs along the width of the beam (i.e., 
chord of the wing) where the shell and beam intersect. Along these lines of 
contact the wing imparts both shear and moment reaction loads to the fuselage. 

If the beam is assumed to pass through the center of the shell, the reactions of 
the beam on the shell are an in-plane load and a moment. In this fully coupled 
structural model the effects of stiffeners in the form of circumferential 
frames or longitudinal stringers as well as the effects of loading on the shell 
imposed by the internal fluid are neglected. The latter effects are discussed 
further in Chapter IV. These effects are considered separately using a model 
of a frame stiffened shell that is described at the end of this section. 

2. Unstiffened Fuselage 

Figure III. 8a, b, and c shows the geometric idealization for the 
analysis of the wing- fuselage system. The wing, modelled as a beam of length 
L^, passes through the diameter of simply supported cylindrical shell of 
length (w and f subscripts respectively denote wing and fuselage) . 

Stiffening frames and stringers are not included in the analytical model that 
follows; however, the effects of frames are examined in Section III.C. As 
indicated in the wing free body, the structural effect of the shell is modelled 
as a set of unknown normal- to -be am forces and moment reactions at the two 
points of contact. For the shell analysis these same forces and moments are 
assumed smeared over, or distributed along, the width of the beam (root chord 
in an actual aircraft) . 

Displacements and rotations at the two contact locations may be 
independently expressed for beam and shell in terms of respective structural 
admittance matrices and the unknown load vector (for the beam the known 
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exciting force F also takes part) . Finally equating displacements and 
o 

rotations for the two systems at their two points of intersection one obtains 

a set of 4 simultaneous equations (2 degrees of freedom x 2 points of contact) 

from which the unknown forces may be calculated. For the shell we have chosen 

the center chord point z in the displacement calculation for continuity with 

o 

the beam motion. Thus the symmetric coupled system becomes. 


c n Kc n 


C 12- KC 12 


C 12- KC 13 


C 14 -Kc 14 


c„ -Kc„_ 
22 22 


C 23 +KC 14 


C 24 +KC 13 


C 33 -KC 33 


C 34 -KC 34 


Mj/a 


C 24- K0 33 




where the constants c , c , etc., form the admittance matrix of the wing, and 
* * ^ ^ ^ 

c„ , c etc., that of the shell. These are all given below. The nondimensional 
11 12 

parameter K gives some indication of the relative importance of the beam and of 
the shell in the combined admittance matrix. It is given by 


4(l-vJ) 




Here L _ , as all wavenumbers and other spatial quantities except a, h, and I 
f , w 

in what follows, appears normalized by the shell radius a. 1^, E^ stand 

respectively for the moment of inertia of the beam about its cross-sectional 

horizontal axis, and E for the beam's Young's modulus. E_, v_, h, and a are 

w t r 

shell values of Young's modulus, Poisson's ratio, thickness and radius. 


functions 


The beam constants are expressed in terms of free-free beam eigen- 
19 


V* 1 = 


sinhk. x+sink. x+Cfcoshk^. x+cosk. x) 
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where r is the beam's radius of gyration (normalized by a) for its cross- 
9 

sectional horizontal axis. 

The beam system's constants are 
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The shell constants may be obtained from analysis for thin cylindrical 
shells in Ref. 13, p. 237. The shell system’s constants are 
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and where stands for the shell loss factor. 


C. Model of Frame-Stiffened Cylinder 

The external fuselage is a shell that is stiffened by both 
circumferential frames and axial stringers. In the laboratory test cylinder 
discussed in Section V the shell is coupled to a beam through either of these 
stiffeners by bolt attachments. The analytical model of the coupled beam and 
shell discussed in the previous section does not include internal stiffeners. 

In order to examine the effects of these structural elements, the structural 
response field of a frame stiffened cylindrical shell with internal fluid in 
response to localized circumferential force and moment excitations is examined. 
These excitations can be viewed as providing influence functions (i.e. , pressure 
per shell load) for the coupled problem of a shell and wing. 

The analytical model of the framed stiffened cylinder assumes the shell 
to be simply supported. The model allows the interior fluid loading to be 
included in calculating the shell response. The frames are modeled as rings. 
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having three translational degrees of freedom as well as rotation aligned with 
the circumferential direction. The impedance matrix for the frame dynamics is 
given in Ref. 20. Coupling between the shell and frames is assumed to occur 
through line reaction loads at each frame location. Excitation loads are 
applied at the centroid of the cross section of the driven frame. Under these 
assumptions each circumferential mode of the shell couples to the frames and 
the reaction loads for a given circumferential mode are obtained by solving a 
system of equations equal to 4 (frame degrees of freedom) times the number of 
frames. The terms in each equation involve a sum over the axial modes of the 
shell. Once the frame reactions are known, they are applied one-by-one to the 
shell to calculate the response. The contribution of any particular frame to 
the response can therefore be examined by omitting it and noting the change in 
the response. As discussed in Chapter V, this procedure is useful for examining 
the importance of specific structural elements to the response and interior 
sound field. 
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Fig. III.l - Schematic of Wing Structural Components 
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Fig. III. 3a - Mathematical model of skin-air/fuel layer interaction. 
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Fig. III. 3b - Mathematical model for skin-rib interaction. 
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y = k fy 



Fi-9* III. 7 - Mathematical model of skin-spar interaction. 
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(b): cylinder free body (c)wing free body 


Fig. III.8a,b,c - Wing- fuselage structural interaction model. 
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IV. ACOUSTIC MODELS FOR CABIN INTERIOR 


A. The Sound Field in a Fuselage Modelled as a Cylindrical Shell 
Terminated by Bulkheads 

1. Introduction 

The sound field in an aircraft fuselage, modelled as a cylindrical 
shell vibrating in response to structureborne noise incident along its wings, 
is formulated for two types of acoustical treatment. In both cases the 
cylindrical shell is assumed to be simply supported and terminated by bulk- 
heads, the wings being located in the horizontal plane of symmetry. Both 
analyses assume that the internal source impedance of the shell modes is so 
large compared to the impedance of the corresponding room modes, that they 
(the shell modes) act effectively as constant-velocity sound sources. Under 
these circumstances, pressure peaks are associated with antiresonance of the 
air volume, as well of course as with shell resonances. 

First, in Section B the sound field inside a cylindrical shell with zero 
surface absorption is calculated as the sum of two pressure components, one 
describing the contribution of the cylindrical part of the boundary, and the 
other that due to the presence of the two end bulkheads. In Section C.l 
this solution is generalized in approximate form to account for small surface 
absorption. In the theory each of the pressure components contains parameters 
through which bulkhead and radial absorption may be independently controlled. 

The result derived here, however, takes into account only the latter effect. 
Inclusion of the former axial effect would require a parallel though somewhat 
simpler analysis which would complicate final expressions needlessly. In 
Section C.2 absorption parameters used in B.l are related to fuselage rever- 
beration time. Finally, Section D is a discussion of the theoretical modal 
density for both shell structural response and internal acoustic field. 

B. The Sound Field in a Vibrating Cylindrical Shell with no Sound 
Absorption 

1. Formulation of the Problem 

The radial shell response is in the form of a double, Fourier 

series. 


IV- 1 


w(4>,z) = exp(-iu)t) E w sinn<J> sin — 

m,n mn L 


(IV-1) 


Here z is measured from the one shell extremity -and the circumferential angle 
♦ is from the plane containing the wings. In the more general case where the 
wings are not restricted to lie in a plane of symmetry, E,.. IV-1 would have to 
be supplemented with a cosntj) series. ' 

The sound pressure satisfies the steady-state wave equation in cylindrical 
coordinates 


(V 2 +k 2 )p = 0 


and the boundary conditions 
3p 
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A solution satisfying the former boundary condition is readily constructed, 
by inspection, but being matched to Eg. IV-1, it necessarily displays nodal 
planes coinciding with the bulkheads. This solution taken by itself is 
incompatible with the latter boundary condition in Eq. IV-3. We shall there- 
fore construct the pressure field by combining two pressure components. 

(IV-4) 
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satisfying respectively the boundary conditions 
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Thia procedure can. if desired, be modified to account for a non-zero specified 


bulkhead motion, viz., E u (r)sinn(fc. 
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2 . Pressure Field Radiated by the Vibrating Shell 


Eq. III-l admits solutions in the form of Bessel functions of the 
first and second kind. However, the latter display a singularity on the axis 
r=0. Consequently, in the absence of significant absorption, only functions 
of the first kind are required. Taking w positive outward, it can be verified 
that the solution which satisfies Eq. IV-5 is of the form 
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where the radial wavenumber is 
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k = [k 2 -(m7r/L) 2 ] i 
m 


(IV-8) 


This pressure component describes the sound field in an infinitely 
periodic shell. If the shell is provided with a sound absorptive core, the 
Bessel functions in Eq. IV- 7 are replaced by Hankel functions. Both solutions 
predict vanishing pressures in the planes z=0 and L. Some analyses even though 
intended for finite shells, erroneously omit the pressure component p required 

2i 

to account for the bulkheads, whether the space is anechoic or not. 


The axial pressure gradient in the plane of the bulkheads is 
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The pressure component in Eq. IV- 7 displays resonances at frequencies 

where the product k^a equals a root of the Bessel functions. As explained in 

the introduction, the pressure peaks generated by a high-^internal impedance 

source are associated with antiresonances at the zeros v a of the Bessel 

22 ns 

function derivatives 
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There is therefore a triple infinity of room resonances and anti resonances . 

For k = 0, 
m 

f a = cm 

om “ 2L (IV-11) 

is an antiresonance for n=0. This wave harmonic is excited only if the wing 
plane does not contain the cylindrical axis. 

It is convenient to distinguish between odd and even axial modes by 

selecting a new origin of coordinates in the midplane, at z = L/2. The new 
axial coordinate is 

z' = z - ih 

The pressure in Eq. IV-7 now becomes 
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3 * P ressure Field Generated by the' Bulkheads 

We now construct the pressure component which satisfies Eq. IV-6, 
where dp ± /^ z is given by Eq. IV-9. Noting that odd values of m in Eq. IV-7 
correspond to the even modes, i.e., modes symmetrical about the midplane, 
while even values of m correspond to modes which are antisymmetric about this 
midplane, the pressure component we are seeking is of the form 
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where 
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and where, as in Eq. IV-11, y a is a root of j'. 

ns n 

Consequently, the introduction of p 2 does not alter the boundary condi- 
tion satisfied at r=a by p r Differentiating Eq. IV-9 with respect toz', 
one can satisfy the boundary condition in Eq. IV-6 by expressing each m,n 
term in Eq. IV— 9 to a Fourier-Bessel series: 
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Multiplying both sides of this equation by J (y r) and integrating allows 

n nS 23 24 

one to solve for the coefficients of the Fourier-Bessel series: ' 
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where the integral I is 
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Combining Eqs. IV-13, 16 , and 17, this yields an explicit expression for the 
bulkhead- reflected pressure: 
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m 1 § 3 , • • • 


sink z 
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cos(k L/2) 
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m 2,4,... 


This does not introduce any additional resonances or antiresonances, 


The resultant pressure is finally obtained by adding Eqs. IV-12 and 18 
as indicated in Eq. IV-4. 
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c * Sound Field in a Vibrating Cylindrical Shell with Small Sound 
Absorption 

1 • Generalization of Section B Theory 

We define the cylinder internal pressure p* and associated 
Green's function G as the absorptive counterparts of those in Section B 
(there defined without asterisk though G is not used explicitly) . With non- 
zero surface admittance K (Ref. 1, Eq. 24) their boundary- value problems are, 
respectively , 
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so that 
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In terms of room modes appropriate to the cylindrical geometry the 
solution for G may be constructed for small values of K: 
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where r,r Q are normalized by the' radius a and y* s are the eigenvalues of the 

radially absorptive boundary condition , a modal form of which (Eq. IV-25 below) 

is given shortly. For small K y may be assumed equal to the real zero y 
, , ns ns 

plus a small complex part; viz.. 
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Y - Y +ie (IV-23) 

ns ns 

For e=0 G* as given by Eq. IV-22 collapses to the room mode form of G , 

the Green's function for the nonabsorptive case which contains an infinite anti- 

2 2 2 

resonance whenever the denominator vanishes, i.e., whenever (ka) - (nma/L) -y ” 0, 

ns 

or at frequencies given by, 

£ 2 irc/a 

J I nnra \ 2 2 

/ \-r) 

With absorption, however, the corresponding value of this quantity is 
always finite, since then 
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We now calculate e in terms of K. The r = 1 homogeneous boundary condi- 
tion on G (Eq. IV-20b) implies that 
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Since J (y* ) - icJ (Y ), an d J (y* ) = J (y )+0(e 2 ) it follows that 
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so that the right side of Eq. IV-24 becomes 
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The p^ component of pressure with nonzero absorption may now be obtained 

★ 

in terms of a triple sum of room modes by substitution of G with the above 
result for e: 
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or 
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In Section B, appears in a different form, however. There the more 
economical double sum solution in Eq. III-7 was the end. result of an implicit 
construction of the Green's function to satisfy the radial boundary condition 
with no absorption. Letting e=0 in Eq. IV-27, it follows therefore that the 
following identity must hold for any m, n, ka combination 
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One may now approximate the right side of Eq. IV-28b for e « 1. The 
result is the term on the left side of Eq. IV-29 plus the following correction 


+2ier (1-r) I 
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(IV- 30) 


Away from the antiresonances then, where no single term in the s sum 
dominates, the eigenfunction for the double-sum solution from Section B may 
still be used as the right side of Eq. IV-29 after neglecting its 0(e) correction 
given in Eq. IV-30. 

At the antiresonances, where the s=s term for which Eq. IV-24 is satisfied 
dominates the sum, we find that the latter may then be approximated to yield that 
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which upon substitution of e in terms of K (Eq. IV-26) , and use of Bessel's 
"22 2 

equation (J /(n -y ) = J /y ), gives that 
n ns n ns 
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In summary then, the absorptive double- sum model for p^ could be kept as 
given in Section B away from the antiresonances; and near them, where the term 

on the right side of Eq. IV-29 becomes infinite, the latter can be replaced by 
the expression on the right side of Eq. IV-32, i.e., a finite quantity. A 
practical payoff of the approximations is that the absorptive case continues 
then to have Bessel functions of real arguments, thus allowing use of well-known 
efficient numerical techniques for their evaluation. 

* 

So far we have dealt only with the p^ component for pressure, which 

appeared in Section B in terms of a constructed Green* s function and which we 

could not therefore immediately generalize to account for small radial 

absorption. The analysis first had to invoke the room-mode formulation, for 

which the generalization to include absorption is in principle always simple. 

T ^ e P 9 solution given in Eq. IV-18, however, appears already cast in room-mode 

* 

form. The corresponding p solution is thus simply given by that same expression 

* ^ 

with Y replaced by Y_ • Furthermore, for convenience of computation and 
ns ns 

. t # * * 
consistency with the above analysis for p / the imaginary part of y may be 

1 ns 

dropped from all Bessel function arguments and kept only in the algebraic 
denominator terms for modes containing an antiresonance. 

It is interesting to note that the p 2 solution in Eq. IV-18 exhibits a 

larger number of antiresonances than does p . It contains all the antiresonances 
. 2 2 

apparent m p , as indicated by the y -k term in the denominator in Eq. IV-18 
(where wavenumbers do not appear normalized by a as they do here) , but also it 
contains others corresponding to zeros of k ns in Eq. IV-14. However, these new 
antiresonances may be interpreted as given by the same criterion which defines 
those of p^. , so long as the latter may be allowed to include the m=0 possibility. 
The reason for omission of m=0 in earlier discussions is, of course, that m=0 
represents a degenerate mode for p^ , which is constructed in terms of sin nmz/L 
eigenfunctions. 
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2. Estimate of Reverberation Time 


It may be useful to relate the notation used above to the symbols 

28 29 

found in standard texts such as Morse and Pierce, and to express the 

attenuation term in terms of the reverberation time. While Eq. IV-23 is 

associated only with the radial wavenumber, being effectively a large-L or 

small-m approximation, room acoustics theory embodies the resultant wavenumber. 

For ease of relating variables to those of Morse and Pierce, we revert to the 

dimensional forms of y , etc., as used in Section IV. B. 

ns 

In cylindrical coordinates, the complex wavenumber is 
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This related to the imaginary term in Eq. IV-23 as 


n„ = - 2 — , n.. « 1 

N y N N 


In Morse's notation. 


28 


y = w /e / ri M = 2 
N N N to. 




N 


In Pierce's notation. 


29 


y N - k(n) 


kB 


nn 


2 

k (n) 


( IV-34) 


The imaginary component of Eq. IV-33 is related to the characteristic decay time 

T of mode N as 
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Im 



k 



Since x N equals the reverberation time divided by 6 in 10 = 13.8, this can be 
written. as 


*2 7.3x10 -2 

Im Y = ^ m 

n 


where f is in Hz and T in sec. 
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The attenuation constant in Eq. IV-3 can be computed from the attenuation 
A^ over the bulkheads and the attenuation A^ over the cylindrical boundary 
(Ref. 28, p. 409) : 


G 


1 

16V 





(IV-35) 


where V is the volume of the cylinder. 


D. Modal Density and Resonance Bandwidth of the Structural Response and 
Sound Field 

Peak pressures are observed at the natural frequencies of the shell 
and at the eigenfrequencies of the room, Eq. IV-10. It is of interest to 
determine the respective contributions of these two families of resonances for 
the parameters of the laboratory test model. 


1 . Structural Resonance Spacing of Cylindrical Shell 


The shell acts effectively as a membrane below its breathing mode 

natural frequency f = 27ra/c , where c is the plate velocity in aluminum. For 

o p p 

membrane modes, stiffeners are not anticipated to raise the natural frequencies 
significantly, i.e., to decrease the modal density. In the frequency range 
above the breathing mode resonance, where the shell responds predominantly in 
flexure, stiffeners are anticipated not only to raise the natural frequencies, 
but to cause the modal density to decrease with increasing frequency, in contrast 
to the modal density of the unstiffened shell which is frequency independent at 

high frequencies. In the low frequency range, below the ring resonance f = 

2 r 
f (1-v ) (= 3.2 kHz for the 0.50 m (20 in.) shell), we shall therefore approximate 

the modal densities of the stiffened shell by those of the unstiffened shell. 
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For the 0.5 m (20 in.) shell, this becomes 

§= 1.3 x 10 _2 f 1/2 , f * 1.5 kHz 
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3.4 x 10 _4 f 


1.5 kHz < f < 1.7 kHz 


The spacing (dN/df) 1 of the resonance peaks is shown in the first column of 
Table IV-1. 

2. Room Mode Antiresonance Spacing 

The modal density of the characteristic frequencies of the modes 
of a cylindrical room are given by (Ref. 28, p. 400) 
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For the dimensions of the 0.50 m (20 in.) shell, this yields 

-1 
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dN -82 -5 -3 

— = 7.6 x 10 f + 3.1 x 10 f + 0(10 ) , 
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The corresponding frequency spacing is tabulated in the second column of 
Table IV-1. 


3. Bandwidth of Pressure Peaks for Structural Response and for 
Acoustic Field 

The bandwidth between the -3 dB points of a peak is rjf* Even 
selecting a conservatively small structural loss factor of n =0.01, the 
bandwidth is seen to exceed the average resonance spacing of the shell modes 
(third column of Table IV-1) . 

The loss factor associated with a room mode is computed from Eq. IV-35 

2 

with = k. We shall conservatively assume that ( n /Y ns a ) << Assuming that 

all interior boundaries are covered with a blanket whose absorption coefficient 
is as indicated in the fourth column of Table IV-1, one computes room bandwidths 
which, once again, are wider than the average resonance, or antiresonance, spacing. 
Since room modes tend to cluster, there may of course exist noticeable pressure 
peaks in frequency ranges sparsely populated by characteristic frequencies, but 
we do not expect marked peaks at the frequency spacings calculated above. It 
should, of course, be recalled that these results are semi-quantitative in that 
we apply an asymptotic small-boundary admittance theory. 

Finally, since boundary absorption limits pressure buildup at these 
antiresonance frequencies, the modelling of the shell as a large-internal 
impedance sound source seems justified. For this same reason it is anticipated 
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that as long as the model shell boundaries are provided with an acoustic blanket, 
while the shell is devoid of a damping coating, marked pressure peaks will be 
associated with structural resonances rather than with the characteristic 
frequencies of the room modes. 


TABLE IV- 1 

Resonance Characteristics of 0.50 m (20 in) Model Shell 

(All Units in Hz) 


f 

Average Resonance Spacing 

-3 dB Resonance 

Bandwidth 

Shell Modes 

Room Modes 

Shell (n=10~ 2 ) 


Room 

500 

3.4 

29 

5 

26 
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t! 
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* — - 

1000 

2.4 

9.1 

10 

46 

(a=0. 7) 

2000 

1.5 

2.7 

20 

46 

P 

II 

o 
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V. CALCULATED RESULTS FOR INTERIOR NOISE 

A. Introduction 

In this chapter we present results for the interior acoustic field 
calculated using the analytical models of the source, structure, and sound 
field developed in the preceding sections. The calculations have been performed 
for parameters representative of two physical structures. The first is a full 
scale aircraft excited by the propeller wake acting on the skin of the wing. 

The second structure is a laboratory model consisting of a shell coupled to a 

uniform beam. Excitation of this model is a shaker applied near the outer end 
of the beam. 

Because our goal in this study has been to examine a variety of parameters 
relevant to the structureborne noise problem, various models (i.e., assumptions) 
have been used to obtain the results presented in this section. 

B . Full Scale Aircraft 

The dimensions used for the estimate of the interior noise field of a 
small aircraft in response to a propeller wake excitation of the wing are given 
in Table V-l. As indicated the chord of the wing is taken to be 1/4 the length 
of the shell, the wing being centered along the shell. The principal assumptions 
invoked to obtain the interior noise estimate are given as follows: 

1. The propeller wake- induced pressures propagate to the fuselage 
along the skin of the wing. 

2. The structural impedance presented to the wing skin by the fuselage 
is large. 

3. The effects of coupling between the wing and fuselage on the side 
away from the excitation are negligible. 

4. The fuselage can be modelled as a uniform shell. 

5. The interior pressure field is consistent with pressure- release 
terminations (i.e., the p^ component discussed in Chapter IV). 

The first assumption defines the structureborne path being considered, 
and relates to the fact that the skins of the wing and fuselage are contiguous. 
The implication of the second assumption is that the wing reaction loads can be 
calculated assuming rigid constraints at the fuselage interface. The reaction 
loads thus calculated can then be applied directly to the shell. The third 
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assumption follows from the second in that the fuselage blocks propagation along 
the wing and any back reaction of the shell on the low impedance skin is 
negligible. The last two assumptions are invoked for convenience in computation. 

The presence of internal fluid is included in this model in evaluating the shell 
response. 

The pressure at the propeller blade passing frequency (283 Hz) and 
normalized to the maximum pressure on the wing in the propeller wake is shown as 
a function of axial distance from the midspan on Fig. V.l. Circumferential 
variations of the pressure have been averaged over angles of 30°, 60° , and 90° 
from the location of the excitation. The parameter of the curves presented is 
the radial distance from the center of the shell. It is noted that the predicted 
level of the interior pressure is comparable to the maximum pressure applied to the 
wing by the propeller wake (i.e., 0 dB) . This result differs from the relatively 
low radiation from discontinuities along the wing predicted earlier (Section III. A. 5) 
primarily because of the high radiation efficiency of the membrane-like modes in 
the (cylindrical) shell. Fortunately as discussed in Section III. A, this 
structural path can be easily interrupted by ribs or other structural discontinu- 
ities. Another feature of the results in Fig. V.l is the smooth variation of the 
average pressure with both radial and axial locations near the excitation. 

C. Laboratory Model 
1. Introduction 

The schematic in Fig. V.2 shows the test model at NASA Langley 
Research Center, the dimensions being given in Table V-I. The physical model 
consists of a thin cylindrical shell reinforced by five stiffeners in each 
direction (i.e., frames and stringers). The box-beam representing the wing passes 
through the shell just above the horizontal and is driven by a shaker near one 
end. The beam and shell materials are, respectively, magnesium and aluminum, both 
assumed in the theory to have a loss factor of .01. The apparatus is suspended by 
a pair of shock chords, each attached to the shell at the uppermost point of the 
flat end cap. 

As described in Section III the analytical model of the coupled beam 
and shell is a somewhat simplified representation of the physical model. In 
particular neither the effects of the stiffeners nor the loading of the internal 
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fluid are included in this model. Both these effects are examined using 
another model and discussed further along in this section. 

2. Detailed Results for the Coupled Beam and Shell Model 

As previously implied, the calculation of the acoustic field by 

the theories of Chapter III.B and IV. A is really a two-step process. First 

the forces and moments at the two wing- fuselage intersection locations are 

solved for from the combined admittance system (Eq. III-50) . Then, the numerical 

values for the reactions so computed (in terms of the drive amplitude F ) become 

o 

drives for the acoustic model of IV. A. As discussed in Chapter IV, in the 
acoustic model these four reactions appear distributed over the root chord. 

In addition to the final acoustic predictions for the laboratory shell- 
beam model, in fact, in order to best interpret them, some intermediate structural 
results are also presented here. Regarding these, we should point out that for 
the unframed shell model rotational admittances remained consistently high over 
the frequency range investigated (600-2000 Hz) , with the result that the two 
moment reactions are always found to be insignificant relative to the two shear 
forces (typically 20-30 dB down) . Furthermore, since their influence 
coefficients in the sound field are found to be of magnitude comparable to 
those of the shear loads, their acoustic contribution is also negligible. We 
therefore omit from discussion the two moment reactions and their part in the 
coupled system of Eq. Ill— 50 (bottom two rows and rightmost two columns) . 

Figure V.3 shows, vs. frequency, the self admittances c^, Kc*^ for the 
beam and shell, respectively. Relative maxima occur for c at 800, 1200, and 
1700 Hz, corresponding to natural frequencies of the free- free beam. The 
fundamental resonance frequency of the beam is 89 Hz. The greater modal 

density of the shell within this frequency range is evident from both the less 

★ 

patterned nature of Kc^ and the milder oscillations for the same loss factor. 

Of particular interest is the fact that the beam and shell admittances are on 
average of comparable magnitude making for efficient coupling of structureborne 
noise between the two structures. 

Figure V.4 similarly shows the two shear reactions F^ , F (normalized by 
Fq) as a function of frequency. Here the efficient coupling mentioned earlier 
is manifested in the force ratios hovering about unity or 0 dB. Also, the fact 
that F^ and F 2 are of comparable magnitude results from both this close 
coupling and the assumed low loss factor. 
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The structural-acoustical influence function" of the shell, that is 
the interior acoustic pressure level per unit tangential force on the shell is 
plotted in Fib. V.5 for two field locations. The predicted oscillations with 
frequency are great, in excess of 40 dB over the band shown, and since the 
modal density of the acoustical space and the shell exceeds the spacing 
between computed points, they are not connected. Pressure peaks occur at 
frequencies corresponding to resonances of the structure, and, with fluid 
loading on the shell neglected in calculating the shell response, at anti- 
resonance frequencies of the interior acoustic space. This latter source of 
peakiness requires further discussion. 

The interior acoustic space being closed exhibits an infinite set of 
resonance frequencies associated with its "room" modes. Between these 
frequencies exist antiresonance frequencies and ignoring acoustic absorption 
a room mode driven at antiresonance yields an infinite impedance. Under this 
circumstance the modal impedance of the acoustic space exceeds that of the shell. 
If the interaction problem is solved rigorously, that is if the acoustic 
impedance is included in the calculation of the shell response then the anti- 
resonant mode is not excited owing to its infinite impedance. However if 
neglected, the contribution of such a mode to the shell response is finite and 
the infinite acoustic impedance yields the artifact of an infinite acoustic 
pressure level at precisely these antiresonance frequencies. Nevertheless, 

interestingly enough, these predicted peaks do reflect at least qualitatively 
a real phenomenon, one that results from the rigorously solved interaction 
problem in the absence of appreciable absorption and structural damping. • 

As the excitation frequency sweeps through each of these modal acoustic 
antiresonance frequencies the acoustic impedance rises in above and then drops 
below the corresponding modal shell impedance and in so doing changes from 
inertial to stiffness-like or vice versa. Thus not precisely at, but in the 
vicinity of, the antiresonance frequency, the reactive components of the 
acoustic and shell impedances will cancel and a peak pressure will result due 
to what may accurately be called a coupled or system resonance. And in the 
absence of energy dissipation these peaks will be infinite. 

Therefore the only way to reduce this second set of peaks to realistic 
values is to incorporate realistic absorption or damping into the model. 

Although, in this report (Chapter IV) and elsewhere, asymptotic approaches 
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have been developed for this purpose, namely small absorption or black body 
31 

radiation models, more rigorous and computationally efficient solutions are 
required, at a minimum for confirmation of these models. 

An example of the effect of fluid loading on the internal pressure field 
in the absence of absorption is given in Fig. V.6. These results are calculated 
for the unframed laboratory shell excited by a point circumferential load at 
midspan, the pressure being evaluated at the same axial location as the load 
and 2 inches (.05 m) in from the surface. In the frequency range indicated 
there are 21 antiresonance frequencies of the internal fluid. However, with 
a structural loss factor of r\ = 0.01 the results for the fully coupled problem 
(i.e., including fluid loading) are found to vary smoothly near the antiresonance 
frequencies. The results of the calculation performed without including fluid 
loading are however strongly influenced by the proximity of the frequency of 
interest to an antiresonance frequency. Consequently these results are shown 
as points. Differences of as much as 20 dB are found between the results of 
the two calculations shown on Fig. V.6. 

With the above background comments Fig. V.7 is presented comparing the 
calculated acoustic pressure normalized to the input force on the beam with 

measurements. Specifically Fig. V.7 is a plot of acoustic pressure p 

2 OS 
normalized by F /ft as predicted by the nonabsorptive model of Section IV. A. 
o 

The measurement curves corresponds to the nominally "beam-horizontal' 1 position 
(as shown in Fig. IV. lc; beam vertical refers to the latter rotated by 90°) . 

One curve is for the beam attached to the shell ring frames, and the other for 
the beam attached to longitudinal stringers. The axial position z is the 
halfway point between the cylinder ends, and the radial position is r = 8 in. 
or 2 in. in from the shell inner surface. For this r,z combination two 
circumferential angular positions are shown in the predicted results: 45° from 

the horizontal plane measured from the drive-side of the beam, and 90°. 
Measurements are for these same values of r,z, and for cf> = 90°. As expected 
the omission of acoustic absorption in these calculations produces levels that 
are high at specific frequencies as well as on-average. Thus in spite of the 
complexity of the mathematical model which for example includes finite length 
shell and acoustic space effects, and the beam-shell interactions, the salient 
features of the measured data are not predicted well. (The effects of the 
omission of stiffening frames is discussed in Section V.C.5.) The lessons 
to be learned from this exercise have been summarized in Section 6. 
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Although broad band analyses are of limited value given the tonal nature 
of the problem at hand they are generally insightful and in this context a 

simple power flow view of the measurements is presented in the following section 
section. 


>totic model and results 


Detailed (deterministic) mathematical models are necessitated by 
the fact that the major problem associated with prop-fan induced cabin noise 
is tonal in nature. Nevertheless, even for tonal problems broadband or 
statistical analyses are often useful in providing additional guidance and 
insight. For example let us consider the implications of balance of power for 
the set of shell experiments being investigated. 

Power is introduced into the system via the shaker at the termination of 
the beam (wing) . The amount of power is given by 


7T . 
XIX 


- F 2 Re(Y) 


(V-l) 


where F is the rms force level and Re (Y) the real component of the drive point 
admittance of the beam. If we ignore finite length effects on Y, in other 
words if we ignore the effects of reflections from the shell-beam interface 
as well as the far end of the beam, then it can be shown that Re(Y) = k f /u)j 
where k f is the wavenumber of freely propagating flexural waves in the beam and 
its locally reacting inertial impedance. 


Let us suppose that the dominant dissipation mechanism is the absorption 
of acoustic energy in the "blanket" that is located along the inner surface of 
the shell. That is, it is assumed that under steady state conditions the power 
flow is from mechanical power in the beam to acoustic power within the 
cylindrical space (without appreciable losses) and that this power is then 
dissipated along the interior of the shell surface by the blanket. This may 
be expressed as follows 


n in = n out = I a S 


(V-2) 


where I q is the space averaged acoustic intensity at the shell surface and 
S - 2iraL the surface area of the shell. Eq. V-2 may now be used to obtain an 
upper bound (in view of the above assumptions) estimate of the interior pressure 
levels. Taking the plane wave admittance for the acoustic space, I = p 2 /pc, 
where p is the space averaged rms pressure. Therefore, 
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or 


i P/F | 2 = (p/p b ) (k f /k) (A b S) 


(V-3) 


|pa 2 /F| 2 = (2Ti 2 ) -1 (p/p b ) (f^fJ^a/DA^Ah (V-4) 

where p and are the mass densities of the acoustic medium (air) and the 

beam (magnesium) respectively, k — oi/c is the wavenumber of acoustic waves. 

Ah and A g are the cross-sectional areas of the beam and shell, and f is the 

critical frequency for flexural vibrations in the beam (the frequency at which 

the wavelength of freely propagating flexural waves in the beam matches the 

acoustic wavelength) . For the shell experiments performed in the laboratory, 

and with the beam "horizontal", f = 238 Hz and we qet 

c 

20 log | P/ r| (re: 1 ft" 2 ) = - 14 - 5 log (f(Hz)] (V - 5 ) 

(Since is proportional to the radius of gyration of the beam, differences 
between horizontal and vertical positioning of the tested beam are expected 
to be less than 3 dB.) Eq. 5 is compared in Fig. V.8 to an envelope of 
experimental data constructed from the six curves corresponding to measurement 
at <J> = 0°, 90°, and 180° , for both ring and stringer shell-beam attachment. 

Also shown is the three point average (<}> = 45° , 90°, and 135°) of the acoustic 
field as predicted using the methods of Chapters III.B and IV. A. The favorable 
comparison between the power balance calculation and the data is gratifying 
and it is noted that the parametric dependence of the prediction (Eq. V-4) 
implies that for this particular model the acoustic field, frequency averaged 
is relatively insensitive to minor structural changes such as beam-shell 
attachment details or beam orientation. 

4. Effect of Structural Damping of the Shell 

As present in the wing model discussed in Chapter III, the 
effect of adding structural damping to the wing is to increase the propagation 
loss and thereby to decrease the excitation of the fuselage. The effects on 
interior noise achieved by structurally damping the fuselage however are 
somewhat less obvious. Shown in Fig. V.9 is the pressure at a point inside 
the unframed test shell for both small (H = 0.01) and large (B = 1.0) damping. 
The excitation is a point circumferential load (shear), applied at a half way 
point between the shell ends. In the case of large damping, the acoustic 
source can be viewed as being the structural near field of the excitation. 
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and the shell away from the drive point acts as a locally-reacting high- 
impedance boundary. For low damping, the shell away from the drive point 
responds as a wavebearing boundary. The results on Fig. V.9 show that 
comparable peak levels are obtained for both values of damping, indicating 
that radiation from the structural nearfield of the excitation is important 
for the unframed shell. This being the case the interior noise level averaged 
over frequency is not expected to be sensitive to the structural damping in 
the fuselage. However as indicated on Fig. V.9 a general statement regarding 
the effect of damping at a specific frequency and location can not be made. 

5. Effect of Shell Ring-Frame Stiffeners 

Results have been calculated for the test cylinder stiffened by 
five frames located along the shell at intervals of eight inches (see Fig. V.10). 

A structural loss factor of 0.01 is assumed. The dimensions of the frames are 
indicated in Fig. V.lOc. Results are calculated every 50 Hz in the middle 
frequency range from 800 to 1200 Hz. For the purpose of examining the effect 
of frames on the interior acoustic field, the shell is assumed to be excited 
by a single localized load on the frame located midway between the simple 
supports. The resulting symmetry permits the response to be calculated using 
only the odd numbered axial modes (i.e., sin imrx/L, m = 1,3,5,...). Additionally, 
the symmetry implies that the reactions are identical for the frames symmetrically 
located about the center frame. Loading on the shell by the interior fluid is 
included in these calculations. 

As a basis for comparison, the pressure at a 'radial distance of 0.05 m 
(2 in) in from the skin is shown on Fig. V.ll for the unframed shell excited 
by a circumferential in-plane load applied at the center. Figure V.lla shows 
the pressure 90 degrees away from the drive at axial locations from one end of 
0.2, 0.4, and 0.6 m (8, 16, and 24 in). Substantial variations (e.g., 20 dB) 
are found within 0.2 m (8 in) axially within this frequency range. The peak 

pressure in this frequency range is approximately -20 dB (re: 0.09 m 

-2 

(1 ft )). The circumferential variation between points located 0.2 m (8 in) 
axially from the drive is shown in Fig. V.llb. Here large differences (e.g., 

20 dB) are found at frequencies in the range between points located 45 degrees 
apart. 

Results for the pressure in the framed shell excited by a circumferential 
load are shown on Fig. V.12. Although substantial differences at specific 


j, 
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frequencies are found between these results and those of the unframed shell 
shown on the previous figure , the general frequency dependence of both results 
has an approximate correspondence. Furthermore the peak pressure level 
calculated with frames (i.e., -20 dB) is within 5 dB of that evaluated in the 
absence of frames. Finally, it is observed that the pressure field in the 
stiffened shell varies somewhat more smoothly with circumferential distance 
than the field in the unframed shell. The most likely explanation for the 
smoothness is that the response of the framed shell is dominated by the lower 
order circumferential modes with their corresponding longer structural 
wavelength. 

Figure V.13 shows the pressure field in response to a moment excitation 
of 1.4 N-m (1 ft-lb) . Although the respective levels are shifted the results 
are similar to those obtained with the in-plane loading shown in Fig. V.12. 

As previously described the shell is excited in the analytical model by 
a PPlying the reaction loads of the passive frames as well as of the driven 
frame. In a highly damped structure the interior pressure field would be 
expected to be dominated by the loads applied by the driven frame. In an 
undamped structure however the driven frame reactions would not be expected 
to dominate the response. The importance of the driven frame for the lightly 
damped test cylinder is indicated in Fig. V.14. Here the pressure at 0.4 m 
16 in) in response to an in-plane load is shown both with and without the 
reactions of the driven frame included. There is no clear trend evident in 
these results. At some frequencies (e.g. , 1000 Hz), the pressure is much 
different without the driven frame. The conclusion at these frequencies is 
that the result is dominated by the pressure induced by the drive frame loads. 
At other frequencies however, the two results are not much different suggesting 
that the contributions of the driven and passive frames are comparable. 

6. Conclusions 

Based on the results presented in this Chapter the following 
conclusions are drawn: 

a. In the absence of significant scattering from wing ribs or 
other major structural discontinuities, the propagation loss of structureborne 
noise along the wing at the propeller BPF is negligible and cabin noise levels 
are comparable to the source pressure levels for a cylindrical fuselage. 
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b. Consistent modelling of the interior fluid requires either 
the inclusion of fluid loading in evaluating the shell response or the inclusion 
of boundary absorption in the interior space. 

c. The interior acoustic levels on average are not appreciably 
altered by the presence of structural damping in the shell. 

d. Compared with the pressure field in a uniform shell, the 
effect of frame stiffeners alter the details of the response in frequency and 
space. The average response and peak levels however are minimally changed. 

e. The asymptotic model discussed in Section V.B.3 gives a 
conservative prediction (i.e., somewhat high) for the average pressure level 
measured in the laboratory model. 
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TABLE V-l 


Dimensions and Properties of Aircraft and Laboratory Model 


Fuselage (Shell) 
Dimensions 
Radius (a) 

Length (L^) 

Skin Thickness (h) 
Material 

Mass Density (p f ) 

Modulus of 
Elasticity (E f ) 

Poisson's Ratio (v^) 

Loss Factor 

Wing (Beam) 

Dimensions 


Small ^ 
Aircraft 

1.11 m (43.8 in.) 
17.3 m (692 in.) 
1.7xl(f 3 m (.07 in.) 

2.7xl0 3 kg/m 3 
(2.5xl0 -4 lb-s 2 /in 4 ) 
10 7 

6.9x10 Pa (10 psi) 

.33 

0.01 


Effective Cross-Sectional 

Radius of Gyration (r ) 

9 

Skin 
Spar 


Length (L ) 
(Span) W 


5x10 4 m ( . 02 in. ) 
5x10 2 m (2 in. ) 

2.2 m (86.4 in.)* 


} 


Width (AL ) 
(Chord) 


4.3 m (173 in.) 


Material 


Mass Density (p w ) 


Modulus of 
Elasticity (E^) 


2.7xl0 3 kg/m 3 
(2.5xl0~ 4 lb-s/in 4 ) 

10 7 

6.9x10 Pa (10 psi) 


Loss Factor 


0.05 


Laboratory 

Model 

0.25 m (10 in.) 

1.2 m (48 in.) 

4 . OxlO~ 4 m (.016 in.) 

2.7xl0 3 kg/m 3 

(2. 5xl0~ 4 lb-s 2 /in 4 ) 

.10 7 

6.9x10 Pa (10 psi) 

.33 

0.01 


2.9xl0~ 2 m (1.13 in.) 

1.7 m (69 in. ) 

7.6x10 3 m (3.1 in.) 

1.7xl0 3 kg/m 3 
(1.6xl0~ 4 lb-s/in 4 ) 

4. 5x10^ Pa (6.5x10^ psi) 

0.01 


♦Assumed length from propeller hub to wing-fuselage interface. 
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AXIAL DISTANCE FROM MIDSPAN (ft) 


Fig. V. 1 - Average fuselage interior pressure in response to wing excitation loads as a function of axial 
distance from midspan at 3 radial locations (1 ft = 12 in = .305 m) . 








ST- 



Fig. V. 4 - Beam- shell reaction forces vs. frequency for the laboratory model. 






Z.T-. 



without 


Fig. V.6 - Effects of fluid loading on the 
pressure in the unframed laboratory 
shell excited at midspan by a point 
circumferential load. 




0 


1111 Measurements' envelope 


Theory 3-point average 
Power equation average 


Fig. V. 8 - Three-point average (* = 45°, 90% 135% wing determinstic acoustic model of Chapter IV, 

n in = n out statistica l average, envelope of experimental data of laboratory model at $ = 45°, 

90°, 135° for both stringer and ring-frame attachment. 

1 1 1 I I i l i i 
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800 1000 1200 
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Fig. V. 9 - Pressure level in response to a point circumferential load on the 
unf rained laboratory shell for loss factors of 0.01 (———•) and 

1.0 (x x) . Pressure location is x = 16 in, r = 8 in, and 

<t> = 90° (from drive) (1 in = .025 m) . 
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(c) 

Fig. V.10 - Model of the frame stiffened laboratory cylinder showing 

(a) dimensions, (b) loading, and (c) frame details (1 in = 
. 025 m) . 
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Fig. V. 11 - Pressure at a radial distance of 8 in. in the 
unframed laboratory cylinder in response to a circumfer- 
ential load located midway along the shell, (a) Dependence 
on axial distance from one end at <j) = 90°. (b) Dependence 

on angle from load at x = 16° (1 in = .025 m) . 



FREQUENCY (Hz) 



P/Fa dB (re: 1ft' 



Fig. V. 12 - Pressure at a radial distance of 8 in. in the framed 
>- -10| — laboratory cylinder in response to a circumferential load located 

midway along the shell. (a) Dependence on axial distance from one 
end at <j> = 90°. (b) Dependence on angle from load at x = 16 in 

(1 in = . 025 m) . 
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FREQUENCY (Hz) 


Fig. V.14 - Pressure at r = 8 in, x = -6 in, and <p = 90° for the framed 

laboratory cylinder excited by a circumferential load calculated 
with all frame reactions included and by omitting the reactions 
of the driven frame (1 in = .025 m) . 
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